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Quarterly  Progress  Report  No.  i 
on 

Project  N5orl-07856 
SURVEY 

In  the  advanced  course  which  I have  been  giving  this  term,  I have  been  treating 
the  energy  band  theory.  I had  hoped  to  get  further;  but  it  has  taken  the  entire  term  to 
cover  the  material.  This  includes  just  Chapter  1 of  the  proposed  report  on  "Electronic 
Structure  of  Solids".  Under  the  circumstances,  it  has  seemed  best  to  include  just  this 
chapter  on  energy  bands  in  the  Technical  Report  No.  4,  and  this  report,  including  also 
a bibliography  of  papers  on  energy  bands,  is  being  prepared  for  the  printer,  and  should 
be  sent  out  not  long  after  this  Progress  Report.  I had  hoped  also  to  treat  the  use  of 
Wannier  functions  in  discussing  lattice  imperfections  and  excitons,  and  various  applica- 
tions of  configuration  interaction  in  solids,  including  ferromagnetism.  Since  this  ma- 
terial was  not  covered,  I shall  propose  to  take  it  up  in  a continuation  of  the  same  course 
in  the  fall,  and  shall  write  the  proposed  additional  chapters  of  the  report  on  "Electronic 
Structure  of  Solids"  at  that  time,  bringing  them  out  as  a further  Technical  Report  No.  5 
at  a later  date. 

This  course  on  the  theory  of  energy  bands  has  furnished  an  occasion  for  exam- 
ining the  relative  merits  of  the  various  methods  in  current  use  for  studying  energy 
bands.  The  conclusion  to  which  we  have  been  led  is  that  the  orthogonalized  plane  wave 
method,  perhaps  modified  to  Herman's  "augmented  orthogonalized  plane  wave  method", 
is  the  most  promising  method  at  present  for  accurate  numerical  work.  Parmenter's 
attempt  to  apply  the  tight  binding  approximation  accurately  to  chromium  has  shown  the 
almost  impossible  complexity  to  which  that  method  leads  when  we  try  to  apply  it  accu- 
rately; the  orthogonalized  plane  wave  method  seems  clearly  simpler.  Parmenter  ac- 
cordingly is  looking  into  ways  of  applying  this  method  to  chromium. 

The  great  drawback  of  the  orthogonalized  plane  wave  method  is  that  it  is  not 
readily  adaptable  to  values  of  the  propagation  constant  other  than  a few  special  points 
in  the  unit  cell,  such  as  the  center  of  the  Brlllouin  zone  and  a few  other  points.  The 
reason  is  that  at  these  symmetry  points,  a number  of  plane  waves  must  have  identical 
coefficients,  so  that  the  number  of  distinct  coefficients  is  relatively  small,  even  for  a 
sum  of  a large  number  of  plane  waves,  and  we  do  not  face  a very  difficult  secular  equa- 
tion. A method  of  interpolation  between  these  symmetry  points  is  highly  desirable. 

With  this  in  mind,  it  occurred  to  me  that  the  tight  binding  approximation  might  furnish 
a good  Interpolation  method,  if  the  various  integrals  which  occur  in  it  were  regarded  as 
disposable  parameters,  rather  than  quantities  to  be  calculated  from  wave  functions  and 
potentials.  Dr.  Koster  and  I have  been  examining  this  possibility,  and  it  seems  very 
encouraging.  In  particular,  we  have  carried  through  a study  of  the  diamond  lattice, 
fitting  the  points  determined  by  Herman  by  the  orthogonalized  plane  wave  method,  and 
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our  solution,  which  describes  both  the  valence  and  conduction  band,  seems  very  likely 
to  be  rather  close  to  the  truth.  We  find  that  we  can  get  this  fit  by  using  only  nearest 
and  second  nearest  neighbor  interactions,  and  including  no  three-center  or  overlap  in- 
tegrals; there  are  enough  remaining  integrals  to  use  as  disposable  parameters.  For- 
tunately, the  integrals  derived  by  fitting  Herman's  values  have  entirely  reasonable 
values,  having  the  correct  signs  and  relative  magnitudes,  so  that  this  procedure  seems 
entirely  reasonable.  We  are  carrying  this  method  further,  considering  somewhat  more 
complicated  crystals.  This  work  Is  report  partly  in  the  present  Progress  Report,  part- 
ly in  Technical  Report  No.  4.  We  are  indebted  to  Dr.  F.  Herman,  of  R.  C.A.,  for 
valuable  help  in  this  problem;  he  has  also  been  independently  considering  the  tight  bind- 
ing approximation  for  interpolation  purposes. 

The  work  on  a two-electron  model  of  ferromagnetism,  on  which  Drs.  Statz  and 
Koster  and  I have  reported  in  the  two  preceding  Progress  Reports,  has  been  written  up 
for  submission  to  the  Physical  Review.  Drs.  Statz  and  Koster  give  in  this  Progress 
Report  a two-dimensional  model,  coming  a little  closer  to  the  truth  than  the  one-dimen- 
sional model  which  I presented  in  the  preceding  Progress  Report,  and  which  is  partly 
incorporated  in  the  Physical  Review  paper.  As  for  future  work  on  this  theory,  we  feel 
that  we  want  to  get  closer  to  the  actual  structure  of  a ferromagnetic  metal,  such  as 
nickel,  and  this  demands  knowing  something  about  the  Wannier  functions  and  energy 
bands.  We  a»"e  delaying  work  in  this  direction,  on  account  of  the  feeling  that  the  tight 
binding  interpolation  method,  which  I have  just  been  describing,  may  well  be  adapted 
to  such  energy  bands,  and  to  the  Wannier  functions  to  be  formed  from  them.  Dr.  Ki- 
kuchi,  meanwhile,  is  looking  at  a very  simplified  two -atom  model,  in  the  hope  that  he 
will  find  a ferromagnetic  ground  state  which  may  well  have  analogies  to  the  ground  state 
in  more  complicated  systems. 

Dr.  Pratt  continues  to  work  on  the  problem  of  antiferromagnetism,  which  is 
one  of  the  most  intricate  problems  which  we  face.  His  two  communications,  and  one  of 
Dr.  Koster  on  an  extension  of  Hund's  rule,  all  have  a bearing  on  this  problem,  whose 
status  is  something  of  the  following  sort.  We  have  examined  a number  of  models;  and 
the  difficulty  which  arises  is  that  there  is  a great  tendency  for  a ferromagnetic,  rather 
than  an  antiferromagnetic,  state  to  lie  lower.  This  is  probably  a result  of  the  method 
of  approximation,  but  It  illustrates  the  extreme  sensitivity  of  the  problem,  and  the  need 
of  examining  it  very  carefully  before  we  really  understand  the  validity  of  the  various 
approximations.  In  particular,  Koster  shows  that  Hund's  rule,  stating  that  the  ferro- 
magnetic state  will  lie  lower,  holds  for  any  case,  molecular  as  well  as  atomic,  in 
which  our  configuration  interaction  problem  is  an  ordinary  spin  degeneracy  problem 
between  orthogonal  orbitals.  Probably  as  a result  of  this,  Pratt  finds  that  the  model 
which  he  has  set  up  of  a triatomic  molecule  (MnOMn)  has  a ferromagnetic  state  as 
its  ground  state,  though  we  had  hoped  that  we  should  find  an  antiferromagnetic  ground 
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state,  and  so  an  explanation  of  superexchange.  We  are  not  satisfied,  however,  that 
this  result  is  correct,  for  very  likely  a much  more  extensive  configuration  interaction 
is  needed  to  lead  to  a possible  antiferromagnetic  ground  state. 

The  writer  suspects  that  the  actual  situation  is  the  following.  He  has  pointed 
out  on  several  occasions  that  the  effective  potential  energy  in  which  an  electron  of  + 
spin  moves  may  be  different  from  that  for  an  electron  of  - spin:  on  account  of  exchange 
effects,  the  electron  of  + spin  has  a lower  potential  energy  when  it  is  in  a distribution 
of  charge  of  other  + electrons,  and  a higher  potential  energy  when  surrounded  by  elec- 
trons of  - spin.  If  we  then  have  an  antiferromagnetic  configuration  of  (MnOMn)+>,  the 
spins  on  one  of  the  Mn  ions  (which  will  certainly  all  be  parallel  to  each  other  in  the 
ground  state,  by  Hund's  rule  for  the  Mn  ion)  will  point  in  one  direction,  say  +,  while 
those  on  the  other  Mn  ion  will  point  ir.  the  other  direction,  say  -.  An  electron  of  + 
spin,  then,  will  find  a lower  potential  well  in  one  Mn  ion  than  in  the  other.  This  will 
hold  for  the  electrons  of  oxygen,  as  well  as  for  those  of  the  manganese.  Thus  the  elec- 
trons of  + spin  in  the  oxygen  will  find  themselves  effectively  in  a strong  field,  pulling 
them  toward  the  manganese  ion  with  ♦ spin,  while  the  electrons  of  - spin  in  the  oxygen 
will  find  themselves  in  an  equal  field  in  the  opposite  direction.  These  electrons  will  be 
polarized  by  these  fields,  just  as  if  the  fields  were  ordinary  electrostatic  fields,  and 
this  will  reduce  the  energy  of  the  molecule.  A corresponding  polarization  will  not  occur 
in  the  ferromagnetic  case,  for  then  the  electrons  of  each  spin,  in  the  oxygen,  will  find 
themselves  effectively  in  zero  field.  The  writer  suspects  that  it  is  the  lowering  energy 
on  accobnt  of  this  polarization  which  stabilizes  the  antiferromagnetic  state. 

This  picture  of  the  role  of  the  oxygen  atoms  supplements  the  remarks  made 
by  the  writer  several  years  ago,  regarding  the  alternating  potential  in  which  an  elec- 
tron finds  itself,  in  an  antiferromagnetic  crystal.  It  was  pointed  out  that  this  alternat- 
ing potential  would  split  the  energy  bands  apart,  the  unit  cell  in  the  crystal  being  twice 
as  great  as  ordinarily  considered,  and  that  in  some  cases  we  could  well  have  the  lower 
band  filled,  the  upper  one  empty,  thus  explaining  the  non-conducting  properties  of  cer- 
tain antiferromagnetic  crystals  which  otherwise  would  be  thought  to  be  conductors. 

But  this  argument  by  itself  does  not  show  why  the  antiferromagnetic  crystal  has  a lower 
energy  than  the  ferromagnetic  arrangement.  Dr.  Pratt  reported  some  time  ago  on  hav- 
ing looked  into  a one -dimensional  model  of  such  a system,  and  he  found  the  ferromag- 
netic state  to  have  the  lower  energy.  It  now  seems  likely  that  the  presence  of  the  oxy- 
gen is  required  to  stabilize  the  antiferromagnetic  state,  in  the  way  described  above, 
though  the  splitting  of  the  energy  bands  would  take  place  just  as  described  earlier  by 
the  writer. 

The  problem  is  then  simply  to  find  the  correct  way  to  describe  this  polarize 
tion  of  the  oxygens.  This  could  be  done  by  a configuration  interaction  with  excited 
states  of  the  oxygen  atom,  but  this  is  a notoriously  poor  way  to  describe  polarization, 
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from  the  point  of  view  of  convergence.  It  is  much  better  to  modify  the  wave  function, 
in  the  way  which  Mr.  Allen  is  investigating,  using  a single  determinant  made  of  orbit- 
als which  are  distorted  as  they  would  be  in  the  presence  of  a field.  In  our  case,  the 
oxygen  orbitals  corresponding  to  electrons  of  one  spin  would  be  distorted  in  one  direc- 
tion, those  corresponding  to  the  other  spin  in  the  other  direction,  so  that  they  would  no 
longer  be  identical  to  each  other.  All  this  is  easier  to  describe  in  simple  language, 
than  it  is  to  carry  out  in  detail.  For  one  thing,  we  are  dealing  with  a singlet  state  of 
the  molecule,  and  we  cannot  say  in  any  proper  way  that  the  spin  of  one  manganese  ion 
is  pointed  in  one  direction,  that  of  the  other  ion  in  the  other  direction;  the  true  wave 
function  is  a combination  of  wave  functions  of  the  two  manganese  ions,  with  a variety 
of  orientations.  The  process  of  building  up  the  wave  function  out  of  determinants  is 
extremely  complicated,  on  account  of  the  many  electrons  involved;  it  is  this  part  of 
the  problem  that  Dr.  Pratt  is  handling  by  his  operator  method.  Finally,  if  the  sus- 
picions regarding  polarization  are  correct,  it  may  be  necessary  to  use  non-orthogonal 
orbitals,  and  the  technique  for  doing  this  properly,  without  neglect  of  the  overlap  in- 
tegrals, which  may  be  the  essential  feature  in  the  problen.,  is  very  complicated.  The 
problems  which  Dr.  Pratt  is  facing  are  not  so  much  problems  of  visualizing  the  quali- 
tative nature  of  the  solution,  as  those  of  formulating  the  solution  in  mathematical  lan- 
guage in  which  we  are  sure  that  we  are  not  throwing  away  the  essential  terms  by  mis- 
take. A simplified  model  which  he  has  been  considering,  since  treating  the  cases 
described  in  this  Progress  Report,  gives  hope  for  thinking  that  we  are  on  the  right 
track  in  our  thinking  about  the  polarization  of  the  oxygen. 

The  various  investigations  into  molecular  structure  are  progressing.  An  im- 
portant result  of  our  consideration  of  configuration  interaction  has  been  the  conviction 
that  there  are  very  few  cases  simple  enough  to  handle  by  a complete  configuration  in- 
teraction. If  we  use  only  a partial  configuration  interaction,  it  then  is  very  important 
to  have  good  molecular  orbitals  to  start  with.  There  are  very  few  cases  in  the  litera- 
ture in  which  really  good  calculations  of  molecular  orbitals,  by  a self-consistent 
method,  have  been  carried  out.  Accordingly,  Dr.  Meckler  has  considered  the  possi- 
bility of  mechanizing  the  problem  of  determining  self-consistent  LCAO  molecular  or- 
bitals, and  has  found  a way  to  do  it  on  the  Whirlwind  Digital  Computer,  which  he  de- 
scribed in  the  present  Progress  Report.  If  this  procedure  can  be  actually  put  into 
operation,  it  might  mean  an  important  advance  in  the  problem  of  finding  molecular 
orbitals. 

The  mechanization  of  the  calculation  of  the  integrals  needed  for  the  water 
problem,  by  Koster  and  Schweinler,  is  proceeding.  This  calculation,  Meckler's  on 
ethane,  Kaplan's  on  ammonia,  Barrett's  on  fluorine,  and  Corbato's  on  HF,  all  tie  in 
together  in  the  calculations  of  integrals,  and  general  procedure,  and  are  all  progress- 
ing, though  slowly.  Unfortunately  the  calculation  of  HF  by  Mr.  Merrifield,  reported 
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on  In  preceding  Progress  Reports,  had  to  be  discontinued  at  a point  where  it  had  not 
yielded  ureful  numerical  results;  a continuation  of  this  progi&m  is  being  studied  by 
Mr.  Corbato.  Various  other  somewhat  disconnected  pieces  of  work  are  under  way, 
as  will  be  evident  from  the  Progress  Report. 

There  has  been  one  change  in  the  personnel  of  the  group  since  the  preceding 
Progress  Report:  Dr.  Statz  has  left  the  group  to  join  the  staff  of  the  Raytheon  Manu- 
facturing Company.  Fortunately,  he  will  still  be  in  the  neighborhood,  so  that  he  will 
be  able  to  join  in  the  general  scientific  activities  of  the  group.  Various  other  changes 
in  personnel  are  in  prospect  for  the  summer  and  the  next  academic  year. 


J.  C.  Slater 


1.  A SIMPLIFIED  TIGHT  BINDING  MEi  ttGD  FOR  ENERGY  BANDS 


Among  all  the  methods  which  have  been  used  for  studying  energy  bands,  the 
tight  binding  approximation,  or  Bloch  method,  has  one  outstanding  advantage:  it  gives 
simple  analytical  formulas  for  the  energy  as  a function  of  the  propagation  constant.  It 
has  compensating  disadvantages,  however:  to  carry  it  through  really  rigorously,  tak- 
ing into  account  all  three -center  integrals,  overlap  integrals,  and  other  complicating 
features,  is  practically  impossible,  as  Dr.  Parmenter's  attempt  to  apply  it  to  the  en- 
ergy bands  of  chromium  has  shown.  It  has  occurred  to  the  writers,  however,  that  the 
method  could  be  extremely  useful  as  an  interpolation  method,  for  getting  the  energy  as 
a function  of  propagation  constant  in  a case  where  accurate  values  of  the  energy  were 
known  at  certain  symmetry  points  in  the  reciprocal  space.  This  method  has  been  writ- 
ten up  in  fair  detail  in  the  Technical  Report  No.  4,  which  will  appear  shortly,  and  for 
that  reason  only  a sketch  of  the  work  so  far  done  will  be  given  in  this  Progress  Report. 

Of  all  methods  which  have  been  used  up  to  the  present  for  energy  bands,  the 
orthogonalized  plane  wave  method  seems  to  the  writers  to  be  probably  the  one  best 
adapted  to  getting  accurate  values  of  the  energies.  The  study  which  Herman^  ^ has 
made  of  diamond  and  germanium  gives  one  the  feeling  that  the  results  are  very  reliable; 
particularly  the  indication  that  the  method  is  converging  rapidly  as  more  and  more 
terms  are  added  to  the  series  is  very  encouraging.  On  the  other  hand,  the  orthogonal- 
ized plane  wave  method  is  not  well  adapted  for  calculations  except  at  particular  sym- 
metry points  in  the  reciprocal  space.  The  reason  for  this  is  very  simple:  at  a sym- 
metry point,  many  plane  waves  must  have  identical  coefficients,  so  that  the  number  of 
independent  coefficients  is  far  less  than  the  number  of  plane  waves  which  are  being 
superposed  to  get  an  approximate  solution  of  Schrddinger's  equation.  Thus,  for  in- 
stance, the  greatest  number  of  independent  coefficients  which  Herman  used  in  his  study 
of  diamond  was  16,  so  that  he  had  to  solve  a 16 -by- 16  secular  equation;  but  the  wave 
function  in  this  case  was  made  of  146  plane  waves.  For  a propagation  constant  not 
having  special  symmetry,  we  should  still  need  146  plane  waves  to  get  comparable  ac- 
curacy, but  now  there  would  be  no  special  relations  between  their  coefficients,  and  a 
146-by-146  secular  equation  would  be  required,  which  of  course  would  be  beyond  our 
present  capabilities. 

We  are  then  faced  with  a situation  where  we  have  very  good  energy  values  for 
certain  values  of  the  propagation  constant,  without  an  equally  good  way  to  interpolate 
between  them.  For  Instance,  in  the  cases  of  diamond  and  germanium  mentioned  above, 
we  have  values  for  the  center  of  the  central  Brillouin  zone,  and  for  the  boundaries  of 
the  zone  along  the  100  and  111  directions.  It  has  now  occurred  to  Herman,  ^ as  well 
as  to  the  writers,  that  it  might  be  possible  to  use  a very  much  simplified  tight  binding 
method,  but  to  choose  the  parameters  entering  into  this  method,  not  by  direct  calcula- 
tion from  atomic  wave  functions  and  potentials,  but  by  regarding  them  as  arbitrary 
parameters,  to  be  chosen  so  as  to  fit  the  values  which  were  accurately  known.  We 
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might  then  hope  that  the  solutions  of  the  tight  binding  problem  would  give  acceptable 
interpolation  values  for  the  energy.  Since  we  are  using  the  method  only  for  interpola- 
tion, we  can  afford  to  make  simplifications  which  would  be  entirely  inadmissible  if  we 
were  using  it  for  an  original  calculation:  wc  can  disregard  interaction  integrals  between 
distant  neighbors,  we  can  neglect  all  integrals  except  two-center  integrals,  and  all  over- 
lap integrals.  We  still  keep  enough  terms,  however,  to  have  enough  disposable  con- 
stants to  fit  the  values  which  we  assume  known;  and  of  course  we  must  keep  such  terms 
that  we  have  a good  hope  that  the  calculation  will  agree  with  the  correct  one  to  an  ade- 
quate accuracy. 

To  illustrate  the  method,  we  have  made  a fairly  complete  study  of  the  diamond 
crystal,  fitting  our  constants  to  Herman's  values.  We  start  with  the  two  interpenetrat- 
ing face-centered  cubic  lattices  which  form  the  diamond  crystal.  We  form  Bloch  waves 

from  s,  p , p , p orbitals  on  one  of  the  lattices,  and  similar  waves  from  the  same 

x y z 

orbitals  on  the  other  lattice.  We  have,  then,  eight  Bloch  waves,  and  we  set  up  a secu- 
lar equation  between  these  eight  function  for  a given  value  of  the  propagation  constant. 

In  this  secular  equation,  we  disregard  all  overlap  integrals  between  atomic  orbitals  on 
different  atoms.  The  matrix  components  of  energy  involve  integrals  of  the  form 

H ^(T  - R)  dv.  where  $a  and  ^ are  two  atomic  orbitals,  (for  instance,  an  s 
and  a px  orbital),  R is  one  of  the  displacement  vectors  from  one  atom  to  another,  and 
H is  the  one -electron  Hamiltonian  operator  of  the  periodic  potential  problem.  It  is 
well  known  that  this  integral  can  be  reduced  to  a sum  of  integrals  of  the  form  *1  r) 

V(r  - R')  ^(r  - R)  dv,  where  the  potential  H is  expressed  as  a sum  of  terms  V(r  - R'), 
spherically  symmetrical  potentials  about  the  atoms  located  at  positions  R\  to  an  ap- 
proximation which  is  sufficient  for  the  present  purposes.  The  integral  above  is  a 
three-center  integral,  provided  R and  R'  are  different  from  each  other  and  from  zero. 

We  disregard  all  such  terms,  and  consider  only  the  two-center  integrals  where  R' 
equals  R or  zero.  Then  our  integrals  are  like  those  found  in  a diatomic  molecule.  We 
can  simplify  them  by  expanding  the  $'s  in  terms  of  functions  quantized  with  respect  to 
the  axis  R of  the  pair  of  atoms.  Thus  if  we  are  dealing  with  p orbitals,  we  have  p<r 
and  pir  components  with  respect  to  the  axis.  We  have  non-vanishing  integrals  only  if 
both  +a  and  ^ have  the  same  component  of  angular  momentum  about  the  axis;  that  is, 
if  both  are  a,  or  w,  or  6 functions.  We  have  integrals,  then,  of  types  which  can  be 
symbolized  by  the  azimuthal  quantum  numbers  of  both  orbitals  (denoted  by  s or  p in 
this  case),  and  by  the  component  of  angular  momentum  about  the  axis  (c  or  w in  this 
case).  Such  integrals  can  be  denoted  in  an  obvious  way  as  (ssa),  (spa),  (ppa),  (pp»), 
and  depend  of  course  on  whether  we  are  dealing  with  nearest  neighbors,  next  nearest, 
and  so  on;  this  can  be  denoted  by  subscripts  1,  2,  etc.  It  is  these  various  integrals 
which  we  regard  as  the  disposable  parameters  of  the  problem. 

We  have  set  up  formulas  for  all  matrix  components  of  energy,  diagonal  and 
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non-diagonal,  between  the  various  Bloch  functions  in  terms  of  these  integrals.  We 
then  have,  in  our  case  of  diamond,  an  eight -by -eight  secular  equation  between  the 
Bloch  functions.  This  equation  can  be  easily  solved  analytically  at  the  required  sym- 
metry points,  and  we  can  get  solutions  involving  nothing  worse  than  a biquadratic  equa- 
tion along  the  100  and  11 1 directions.  We  can  now  try  to  fit  all  of  Herman's  values  at 
the  symmetry  points.  We  find  that,  using  nearest  and  second  nearest  neighbors,  we 
do  not  have  quite  enough  disposable  constants  to  fit  all  of  these  values.  However,  we 
can  adjust  the  constants  to  fit  them  all  with  only  small  errors.  In  other  words,  there 
are  relations  between  Herman's  values  which  almost,  though  not  quite,  make  tnem 
compatible  with  a solution  involving  only  second  nearest  neighbors.  When  we  do  this, 
the  resulting  constants  are  Interesting,  in  that  the  integrals  for  second  nearest  neigh- 
bors come  out  much  smaller  than  those  for  nearest  neighbors,  suggesting  that  we 
really  have  a convergent  process.  The  integrals  for  nearest  neighbors  have  the  signs 
and  general  magnitudes  which  we  should  expect;  for  instance,  (ppw),  the  interaction 
between  pv  orbitals  on  the  two  atoms,  comes  out  considerably  smaller  numerically  them 
(pp<r),  the  interaction  of  pa's  on  the  two,  as  we  should  expect  since  the  pa's  overlap 
more. 

When  we  use  these  constants,  we  then  can  get  solutions  for  other  points  of  the 
momentum  space;  we  have  carried  out  solutions  along  the  100  and  111  directions  (ex- 
cept that  we  have  so  far  not  solved  for  the  states  which  require  a biquadratic  equation). 
The  energies  resemble  the  interpolations  which  have  beer,  used  by  Herman  in  his  papers, 
but  differ  in  some  respects;  these  solutions  are  given  in  Technical  Report  No.  4.  The 
qualitative  behavior  of  the  solutions  is  as  given  by  Herman.  Also  the  solutions  in  the 
neighborhood  of  the  degeneracies,  as  in  the  neighborhood  of  the  three-fold  degenerate 
states  at  k = 0,  have  the  behavior  which  we  must  expect  from  symmetry,  and  which 
has  been  described  by  Shockley.  ^ 

The  method  which  we  have’  described  for  diamond  will  be  recognized  as  essen- 
tially equivalent  to  that  proposed  by  Morita.  ^ His  formulas  for  the  matrix  components 
of  energy  are  essentially  like  ours,  except  that  he  does  not  express  his  integrals  in 
terms  of  our  quantities  such  as  (ssa),  etc.  He  finds  his  integrals  from  atomic  orbitals 
and  potentials.  We  have  worked  backward  from  Morita's  published  values  of  his  in- 
tegrals, to  see  the  values  which  he  must  have  found  for  the  various  integrals  which  we 
have  determined  by  fitting  Herman's  values.  In  doing  so,  we  have  observed  what  ap- 
pears to  have  been  an  error  of  some  sort  in  Morita's  integrals.  Specifically,  he  finds 
a numerically  greater  value  for  his  integral  which  we  denote  by  (ppw)  than  for  (ppa), 
contrary  to  our  findings,  and  contrary  to  what  one  would  expect  from  straightforward 
calculation.  This  apparent  error  in  Morita's  calculation  seems  to  explain  another  fea- 
ture of  his  results.  He  found  that  the  lowest  energy  in  the  conduction  band,  for  k = 0, 
came  from  an  s-like  rather  than  from  a p-ilke  state,  contrary  to  Herman's  result.  We, 
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however,  have  been  able  to  find  values  of  the  integrals  which  not  only  put  the  levels  in 
the  same  orJer  as  found  by  Herman,  but  agree  quite  well  quantitatively  with  Herman's 
results.  We  feel,  therefore,  that  the  doubts  about  the  tight  binding  approximation  ex- 
pressed by  Herman  in  the  reference  quoted  are  not  justified. 

A different  sort  of  tight  binding  approach  has  been  carried  out  by  Hall,  ^ 
though  he  does  not  like  to  admit  that  his  method  is  really  a tight  binding  calculation. 

He  starts  with  what  he  calls  equivalent  orbitals.  We  can  describe  these  in  the  following 
fashion.  We  start  with  tetrahedral  directed  orbitals  on  each  of  the  atoms  of  the  crystal, 
directed  toward  the  neighboring  atom.  We  then  make  a symmetric  combination  of  the 
two  directed  orbitals  at  the  two  ends  of  each  bond,  and  call  this  symmetric  combination 
an  equivalent  orbital.  We  have  four  such  equivalent  orbitals  per  unit  cell  of  the  crystal, 
so  that  we  can  make  four  Bloch  combinations  of  them.  We  then  can  set  up  matrix  com- 
ponents of  the  energy  between  these  four  Bloch  functions,  and  can  solve  a four-by-four 
secular  equation  between  them.  This  secular  equation  will  have  four  roots,  which  may 
be  expected  to  describe  the  four  levels  in  the  valence  hand  of  diamond.  It  is  this  secu- 
lar equation  which  Hall  has  considered. 

The  results  which  he  has  found  are  not  in  good  agreement  with  those  of  Her- 
man; in  particular.  Hall,  like  Kimball  and  Hund  and  Mrowka^  with  the  cellular  approxi- 
mation, finds  a doubly  degenerate  state  whose  energy  is  independent  of  the  propagation 
vector,  whereas  Herman  finds  that  the  energy  of  this  state  is  a function  of  k,  and  it  is 
split  in  certain  dJrections;  our  tight  binding  calculation  agrees  entirely  with  Herman  in 
these  respects.  One  way  of  describing  this  shortcoming  of  Hall's  treatment  is  that  he 
takes  into  account  only  interactions  of  nearest  neighbors;  if  he  considered  more  distant 
neighbors,  he  could  get  a suitable  behavior.  In  fact,  Herman^  has  found  it  possible 
to  get  a good  fit  with  his  values  by  making  up  equivalent  orbitals  like  Hall's,  and  using 
sufficiently  distant  neighbor  interactions.  The  interactions  must  be  assumed  significant 
out  to  considerably  more  than  the  second  nearest  neighbors,  however,  which  are  the 
only  ones  which  we  need  use  in  our  method  of  formulation,  in  which  we  get  the  conduc- 
tion band  as  well  as  the  valence  band  from  and  eight-by-eight  secular  equation. 

One  can  easily  convert  Hall's  calculation  into  one  equivalent  to  ours,  and  it  is 
instructive  to  see  how  this  is  to  be  done.  Just  as  Hall  has  set  up  four  Bloch  functions 
from  the  symmetrical  equivalent  orbitals,  we  could  maks  antisymmetric  equivalent  or- 
bitals, and  make  four  Bloch  functions  from  them.  The  four-by-four  secular  equation 
resulting  from  these  would  give  a description  of  the  conduction  band.  But  we  should 
find  that  there  were  non-diagonal  matrix  components  of  energy  between  the  four  Bloch 
functions  formed  from  the  symmetric  orbitals,  and  the  four  formed  from  the  antisym- 
metric orbitals.  These  non-diagonal  matrix  components  vanish  for  k = 0,  but  become 
larger  as  we  go  out  in  the  Brillouin  zone.  If  we  now  consider  all  eight  of  these  Bloch 
functions,  and  the  eight -by -eight  secular  equation  between  them,  we  can  show  that  it  is 
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equivalent  to  the  problem  as  we  have  set  it  up.  We  can  get  satisfactory  results  in  this 
case  by  considering  only  nearest  and  second  nearest  neighbors,  while  to  get  satisfactory 
results  from  Hall's  procedure  we  need  terms  coming  from  more  distant  neighbors.  The 
easiest  way  to  remove  the  zero-width  bands,  then,  is  to  consider  interaction  between 
the  valence  and  conduction  band.  The  point  is  that  except  at  k = 0,  the  wave  functions 
of  the  valence  band  are  not  constructed  entirely  from  the  symmetric  combinations  of 
orbitals  on  the  two  lattices,  and  the  conduction  band  entirely  from  the  antisymmetric 
combination;  closer  examination  of  the  wave  functions  shows  that  the  phase  difference 
between  the  waves  on  these  two  lattices  changes  from  zero  or  180°  at  k = 0 (for  the 
symmetric  and  antisymmetric  solutions)  to  90°  at  the  edges  of  the  band,  a situation 
which  can  be  described  by  Hall's  method  only  with  difficulty. 

We  are  examining  the  tight  binding  approximation  further,  and  hope  to  come 
back  to  the  relation  between  it  and  the  Wannier  functions  for  overlapping  bands. 
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2.  ENERGY  BANDS  IN  CHROMIUM 


The  original  plan  of  calculating  electronic  energy  bands  in  chromium,  reported 
in  previous  Progress  Reports,  involved  expanding  the  valence -electronic  wave  functions 
as  arbitrary  linear  combinations  of  orthogonal! zed  plane  waves  plus  3d  Bloch  functions. 
Each  3d  Bloch  function  was  to  be  a so-called  Bloch  sum  of  3d  atomic  orbitals  associated 
with  the  isolated  chromium  atom.  In  calculating  matrix  elements  of  the  crystal  poten- 
tial between  two  3d  Bloch  waves,  there  arise  terms  involving  three-center  atomic  in- 
tegrals (the  one-electron  integral  of  a product  of  an  atomic  potential  centered  on  one 
lattice  site  and  two  3d  atomic  wave  functions  centered  on  two  different  lattice  sites). 
Because  of  the  difficulty  in  accurately  calculating  three-center  integrals,  the  above 
scheme  is  not  feasible  as  it  stands. 

Let  us  suppose,  however,  that  the  3d  atomic  orbitals  are  modified  in  such  a 
way  that  they  vanish  at  distances  greater  than  one -half  the  nearest  neighbor  distance  in 
chromium.  All  three-center  integrals  will  now  vanish,  and  there  will  be  no  difficulty 
in  calculating  matrix  elements  of  the  crystal  potential.  Such  a method  should  be  general- 
ly useful  among  the  transition  elements.  It  has  been  suggested  by  Herman^  that  the 
same  method  (which  he  calls  the  augmented  orthogonalized  plane  wave  method)  should 
be  more  rapidly  convergent  than  the  orthogonalized  plane  wave  method  alone  when  deal- 
ing with  valence  crystals  (such  as  diamond)  or  semiconductors  (such  as  germanium). 

With  regard  to  transition  elements,  the  motivation  of  the  AOPW  method  is  to  represent 
by  orthogonalized  plane  v/aves  that  portion  of  the  crystal  wave  function  resulting  from 
both  4s  atomic  orbitals  and  the  overlapping  portions  of  3d  atomic  orbitals  and  to  repre- 
sent by  Bloch  waves  that  portion  of  the  crystal  wave  function  resulting  from  the  non- 
overlapping portions  of  the  3d  atomic  orbitals. 

The  method  to  be  used  in  modifying  the  3d  atomic  orbitals  is  more  or  less  ar- 
bitrary, of  course,  but  the  following  method,  which  has  been  used  on  chromium,  seems 
physically  reasonable  and  probably  can  be  used  whenever  the  AOPW  method  is  applica- 
ble. Let  us  construct  an  atomic  orbital  by  numerically  integrating  Schr&dinger's  equa- 
tion with  a Hamiltonian  containing  the  crystal  potential,  the  numerical  integration  start- 
ing at  any  nucleus  of  the  crystal.  Since  we  desire  an  orbital  of  definite  angular  momen- 
tum (i  = 2),  we  must  use  the  spherical  average  of  the  crystal  potential  in  Schrodinger's 
equation.  This  spherical  average  can  be  approximated  by  the  isolated-atom  potential 
V(r)  for  r < R and  by  the  constant  V(R)  for  r > R.  The  radius  R is  chosen  such  that 
V(R)  * <V(r)>,  where  <V(r)>  is  the  average  value  of  V(r),  averaged  over  the  unit 
cell  of  the  crystal.  It  seems  reasonable  to  take  the  orbital  equal  to  the  isolated-atom 
3d  orbital  for  r < R;  i.  e. , we  use  the  energy  of  the  isolated -atom  3d  orbital  as  the 
energy  used  in  numerically  integrating  SchrOdinger's  equation.  For  r > R,  the  numeri- 
cally integrated  orbital  will  deviate  sharply  from  the  isolated -atom  wave  function  since 
the  spherical  average  of  the  crystal  potential  is  lower  than  the  isolated-atom  potential. 
This  will  result  in  a radial  node  in  the  orbital  at  some  value  of  r which  is  approximately 
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equal  to  one -half  of  the  nearest  neighbor  distance  in  the  crystal.  (For  chromium  R = 

1.  55  atomic  units,  the  radial  node  in  the  numerically  integrated  orbital  occurs  at  2.80 
atomic  units,  while  one-half  of  the  nearest  neighbor  distance  is  2.  36  atomic  units. ) By 
cutting  off  the  numerically  integrated  orbital  beyond  the  radial  node,  we  have  formed 
a modified  3d  orbital  suitable  for  our  purposes. 
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3.  A TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM 


In  the  previous  two  Progress  Reports,  J.  C.  Slater,  ^ H.  Statz,  ^ and  G.  F. 
Koster^  have  discussed  a two-electron  example  of  ferromagnetism.  In  this  Progress 
Report  we  apply  the  same  methods  to  a simple  two-dimensional  case.  This  case  fur- 
ther confirms  the  necessity  of  degenerate  bands  for  ferromagnetism  and  also  illustrates 
some  of  the  techniques  which  would  be  useful  in  a more  detailed  discussion  of  ferro- 
magnetism. 

A two-electron  function  in  a crystal  may  be  expanded  in  terms  of  products  of 
Wannier  functions  associated  with  the  various  energy  bands  of  the  crystal. 
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Here  i refers  to  the  band  from  which  the  Wannier  function  a^  arises,  r^  and  r^  are 
the  coordinates  of  the  two  electrons,  and  R^  is  one  of  the  primitive  translations  of  the 
crystal.  We  know,  however,  that  the  symmetry  properties  of  a wave  function  depend 
only  on  the  invariances  of  the  Hamiltonian  and  not  on  the  number  of  electrons  the  wave 
function  describes.  We  expect  therefore  that  all  the  discussion  of  the  irreducible  rep- 
resentations of  the  space  groups  by  Seitz^  and  others^  will  be  applicable  here.  In 
particular,  we  are  able  to  define  a if  yector  which  expresses  the  fact  that  our  two- 
electron  wave  function  goes  into  KP  times  Itself  under  translation  through  it  . 

This  allows  us  to  write  our  U . .(^  , ft  ) in  a simpler  form 
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More  can  be  said  since  we  know  that  our  Wannier  functions  may  have  certain  symmetry 
properties.  Let  us  see  how  this  influences  our  F's.  We  restrict  our  attention  to  K * 0. 
(The  following  arguments  could  be  easily  carried  out  for  other  if  values. ) We  must 
first  see  what  effect  some  operation  a of  our  point  group  has  on  the  wave  function  since 
we  know  that  the  wave  functions  for  K * 0 must  form  irreducible  representations  of  the 
point  group  of  the  crystal.  It  can  easily  be  seen  that 


nj»(r 


1*  r2^  = ^m^  °®i^rl  " *^p)  oa/r2  ^p  + 

i.  J.p.  m 


Here  aa^r)  is  the  Wannier  function  which  results  from  applying  the  operation  a to 
a^(r).  o"  1 is  the  primitive  translation  arrived  at  by  applying  ct“  1 to  ftm.  Since  we 
know  the  transformation  properties  of  our  Wannier  functions  from  Eq.  (3.  3)  we  are 
able  to  get  relations  between  our  F(^m)’s.  This  is  most  easily  seen  by  a simple  example. 
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(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


Assume  that  our  two-electron  wave  function  is  derived  from  the  Wannier  functions  of  a 
single  non -degene rate  band  for  which  oa(i^)  = a(r)  is  true  for  all  o.  Let  us  look  for 
the  two-electron  wave  function  that  forms  a completely  symmetric  representation  of 
our  point  group  (a<|i(r^,  r^)  = r^);  for  all  a).  By  applying  these  relations  to  Eq. 

(3.  3)  we  get 


I F(jtm)a(r  - if )a(r2  - it  ♦ itm)  = I F(o"  1 if  ) a(?  - if - if  + i*m) 

p,  m r r p,  m 


therefore  F(ifm)  = F(a"  1 ft^) 


(3.4) 


m 


The  additional  symmetry  property  of  our  two -electron  wave  function  concerns  its  value 
of  total  spin.  This  is  determined  by  the  relation 

F(-itm)  = t F(ifm)  (3-5) 


The  - sign  refers  to  the  triplet  and  the  + sign  refers  to  the  singlet. 

Let  us  now  apply  the  prece  .ing  discussion  to  a two-dimensional  square  lattice. 
We  wish  to  determine  our  F's  by  minimizing  the  expectation  value  of  a two-electron 
Hamiltonian 


II  = Hj  + H2  + g12  (3.6) 

where  Hj  and  H2  sure  the  one-electron  Hamiltonians  of  our  two  electrons  in  a periodic 
potential  and  gJ2  is  the  Coulomb  repulsion  between  two  electrons  or  two  holes.  We 
study  first  a non-d*»generate  band  for  which  oa(r^)  = a(?)  for  all  a of  the  point  group. 
This  gives  rise  to  the  secular  equation  (Ref.  1,  Eq.  (1.6)) 

I 2*(if J F(Itn  + if J + 6(lfn,  0)  Vq  = EF(ifn)  (3.7) 

P 

Here 

/(itp)  » /aCj)  Hj  a(rt  - itp) 

/*a(r.)  a(r.)  a(r“2)  a(r“2) 

Vn  = 2/ — i = 1-  • (aalaa) 

° J r „ 

12 
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(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


Fig.  3-1 

Square  lattice  in  Rj  - R2  space.  The  figure  gives  the  notation 
for  the  g lattice  points  and  the  various  "planes"  of  reflection. 

In  this  secular  equation  the  assumption  has  been  made  that  two  electrons  interact  only 
when  they  are  at  the  same  lattice  site.  We  also  assume  /-(it  ) 2 0 except  for  nearest 
neighbor  interactions.  For  boundary  conditions  in  relative  coordinate  space  we  make 
the  assumption  that  F(itn)  vanishes  on  the  sides  of  a square  that  contains  nine  lattice 
sites.  Fig.  3-1  contains  a diagram  of  this  region  of  relative  coordinate  space  as  well 
as  labels  for  the  various  lattice  points  under  consideration.  Table  3-1  gives  the  vari- 
ous irreducible  representations  of  group  C4y  (the  point  group  of  a square  lattice).  The 
r's  are  reflections  through  the  "planes"  indicated  in  Fig.  3-1.  C4  is  a clockwise  rota- 
tion through  90°  and  = C4^. 

Table  3-1 

Character  table  for  the  group  of  the  square  lattice  C4y 


\Qpe  ration 
Representatioft\ 

E 

C2 

c4,c4 

°vl ' ®V 

°dl*  °d2 

A1 

1 

1 

1 

1 

1 

A2 

1 

1 

1 

- 1 

- 1 

B1 

1 

1 

- 1 

1 

- 1 

B2 

1 

1 

- 1 

- 1 

1 

E 

2 

- 2 

0 

0 

0 
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(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


We  will  now  discuss  in  detail  one  of  the  symmetries  of  twc-electron  wave 
functions  for  K = 0 under  the  above  assumptions.  Let  us  look  at  those  wave  functions 
with  Aj  symmetry.  Using  the  relation  (3.  4)  we  see  at  once  that 

F(  1,  1)  = F(  1 , - 1)  = F(-  1,-1)  = F(  - 1,1)  = 1 

F(l,  0)  = F(0,  1)  = F(-  1,  0)  = F(0,  - 1)  = | (3.  8) 

F(0,  0)  = a 


Substituting  these  relations  into  secular  equation  (3.  7)  we  find 

a(VQ  - E)  + |(8/)  + 1(0)  = 0 

a(2  / ) + f(-E)  ♦ 1(4/)  = 0 (3.9) 

Z “ 

a(0)  + |(4/)  + 1(-E)  - 0 

therefore 

E3  - V E2  - 32/2E  + 16 / 2V  = 0 (3.10) 

o o ' 


Here  / is  the  nearest  neighbor  interaction  (that  is,  the  matrix  element  of  the  one- 
electron  Hamiltonian  between  a Wannier  function  at  (0,  0)  and  one  at  (0,  1)).  The  equa- 
tion (3.  10)  can  be  solved  for  small  / giving  the  result 


E 

E 


V 

o 


t 4/ 


(3.  11) 


These  states  are  clearly  singlet  states  since  F(^n)  3 + F(-?n).  A similar  discussion 
can  be  carried  out  for  the  other  symmetries.  First  we  write  down  the  relations  be- 
tween the  F's  which  result  from  the  irreducible  representation  which  we  desire  and  the 
symmetry  properties  of  the  Wannier  functions.  We  then  use  these  relations  to  simplify 
the  secular  equation  (3.  7)  and  finally  we  find  approximate  eigenvalues  of  the  secular 
equation  if  these  equations  are  above  ..ie  second  order. 

For  A ^ symmetry  we  find  that  all  the  F's  must  vanish  and  therefore  we  cannot 
find  a state  of  this  symmetrv. 
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(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


For  Bj  symmetry  we  find 


F(l.l)  = F(-l.l)  = F(  1 , - 1)  = F(-  1,  - 1)  = 0 

F(0,  1)  = F(0,  1)  = F(1.0)  = - F(-l,  0)  = £ (3.  12) 

F(0,  0)  = 0 

which  leads  to  the  secular  equation 

|(-E)  = 0;  E = 0 (3.  13) 


This  is  a singlet. 

For  symmetry  we  find 

F(l,  1)  = F(-l.-l)  = - F( - 1,  1 ) = -F(l.-l)  = | 

F(0,  1)  = F(  1 . 0)  = F(-1.0)  = F(0,  - 1)  = 0 (3.14) 

F(0,  0)  = 0 

This  gives  rise  to  the  secular  equation 

|(-E)  = 0;  E = 0 (3.  15) 


once  again  this  is  a singlet. 

For  E symmetry  we  have  a doubly  degenerate  wave  function.  In  this  case  we 
set  up  two  wave  functions 


+y  ■ a(ri  ' Sp>  a(?2  - Sp  * ffm> 


Here  4*  and  4*  transform  as  x and  y do  under  the  operations  of  the  group  of  the 
x y 

square.  Once  again  this  gives  us  relations  between  the  F's 


Fx(l,  1)  = Fx(-  1,  1)  = -Fx{-  1,-1)  = -F^l,  - 1)  = | 

F^l,  1)  = -F3^-  1,1)  = - F^(-  1,-1)  = F^l,  - 1)  - | 


(3.  17) 
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(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


FX(0,  I)  = -F^O,  -I)  = Fy(l,  0)  = -F*(-1.0)  = 

Jl 

FX(1 , 0)  = Fx(-  1,  0)  = Fy(0,  1)  = 1^(0,  - 1)  = 0 (3.  17)  con'd. 

FX(0,  0)  = Fy(0.  0)  = 0 

The  secular  equation  'or  the  x-like  functions  is 


“ (-  E)  ♦ &(4<f)  = 0 

Jl  2 

-^.(2*)  + £(-  E)  = 0 

■Si  2 

with  roots 

E « t lj 2* 


(3.  18) 


Thus  we  get  two  doubly  degenerate  eigenvalues  (a  total  of  four  wave  functions).  These 
states  are  all  triplets  (F(^n)  - - F(-  i$n))*  Fig.  3-2  contains  a graph  of  energy  as  a 


*KI 


Energy  as  a function  of  nearest  neighbor  inter- 
action £ for  the  case  of  non-degenerate  bands. 

function  of  £ for  the  two-electron  states  arising  from  a non-degenerate  band.  We  no- 
tice that  the  singlet  state  of  Aj  symmetry  is  the  lowest  state.  The  fact  that  the  singlet 
lies  lower  than  the  triplet  we  interpret  as  an  absence  of  ferromagnetism. 

We  now  consider  the  case  of  two  degenerate  bands  which  we  describe  in  terms 
of  Wannier  functions  with  x-  and  y-like  symmetry.  ^ We  will  call  them  a(r)  and  b(r ) 
and  their  symmetry  properties  are  given  in  Table  3-2. 
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(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


Table  3-2 

Symmetry  properties  of  Wannier  functions  a and  b 
Operation  b a 


E 

b 

a 

C4 

- a 

b 

C2 

- b 

- a 

C ^ 

^4 

a 

- b 

*dl 

a 

b 

*d2 

- a 

- b 

%1 

b 

- a 

a v2 

- b 

a 

We  write  our  two -electron  function  as 


?2>  ■ £Faa<Rm)  a(?l  ' Bp>  a(?2  • Sp  * Bm> 

* IFbb<Bm>  - Bp>  b<?2  - Bp  * Bm> 

* ZF„b<Bm>  • Bp>  b(?2  • % * Bm> 

* XFb.<B„,>  b<Bl  - Bp>  a<?2 ' Bp  * Km» 


(3.  19) 


If  we  assume  that  there  is  only  interaction  between  electrons  when  they  are  at  the  same 
lattice  site,  these  functions  give  rise  to  the  secular  equations^ 

£J'a(ffp>  Faa(BntBp>  * ‘(Bn’ °>  [Faa(0)(aalaa)  * Fbb<°Xab'ab>]  = EFaa(Bn> 

I 2/b(ffp)  Fbb(^n  + itp)  + 6(itn,  0)  [Fbb(0)(bblbb)  + Faa(0)(ab|ab)j  = EFWBn* 

Z ['a(Bp>  * 'b«M  Fab(Bn*  Bp>  * 6<Bn’  °>  [Fab«»<aa'bb>  * Fba(0)(«bl»b)]  . EFab(S„) 

Z [>a(Sp>  ♦ «b(Sp)]  Fba(Rn  ♦ 5p)  ♦ «Sn,  0)  [Fba(0)(aalbb)  4 Fab(0)(ablab)]  = EFba<Sn) 

P (3.20) 

The  secular  equations  break  into  two  parts,  one  connecting  the  F&a  and  F^b  type  two- 
electron  wave  functions  and  the  other  connecting  the  F b and  Fba  types.  This  splitting 
is  due  to  our  approximation  in  considering  interactions  of  electrons  only  when  they  are 
on  the  same  lattice  site.  In  what  follows  we  make  use  of  these  abbreviations 


-19- 


(A  TWO-ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


V{  = (aalaa)  = (bblbb)  = f alr^)  gJ2  atr^)  a(r2)  di-j  dr2 

J = (ablab)  = Ja^)  b(T2)  g12  alrj)  b^)  d-r  j dT2 

K = (aalbb)  = J a(Tj)  afr^)  gJ2  a^)  bfr^)  dr  j dT£ 

We  also  make  the  assumptions  that  only  those  Fourier  coefficients  of  the  one-electron 

energy  that  represent  nearest  neighbor  interactions  are  non  zero.  We  call  £ the  in- 
teraction between  an  x-like  Wannier  function  and  one  displaced  a unit  translation  in  the 
x direction.  We  call  £ 2 the  interaction  between  an  x-like  Wannier  function  and  one 
displaced  a unit  displacement  in  the  y direction.  A new  feature  comes  into  this  calcu- 
lation of  the  energies  of  the  two-electron  functions  associated  with  the  degenerate  bands. 
In  this  case  it  sometimes  occurs  that  the  symmetry  of  the  wave  function  is  not  sufficient 
to  determine  whether  the  wave  function  is  a singlet  or  a triplet.  When  this  is  the  case 
there  can  be  both  a singlet  and  a triplet  of  the  same  symmetry.  Instead  of  presenting  a 
detailed  discussion  of  the  many  symmetries  and  multiplicites  we  present  the  results  in 
Table  3-3.  The  arguments  that  go  into  deriving  the  tabulated  results  are  exactly  analo- 
gous to  those  already  presented  for  the  non -degenerate  band.  The  first  column  (I)  rep- 
resents the  multiplicity  of  the  state.  The  second  column  (II)  gives  the  symmetry.  The 
third  column  (III)  tells  whether  the  wave  function  comes  from  the  F type  secular 
equations  or  the  F ^ type  secular  equations  (see  Eq.  (3.  20)).  In  the  fourth  column  (IV), 
there  is  listed  the  secular  deterrr  nant  and  in  the  final  column  (V)  is  listed  the  approxi- 
mate eigenvalues  for  £ ^ and  £ 2 small.  In  this  column  we  have  used  the  abbreviation 
that  £^  = c £^. 

For  the  case  of  a two-electron  wave  function  drawn  from  two  degenerate  bands 
it  is  not  possible  to  draw  an  unambiguous  conclusion  as  to  whether  the  triplet  or  the 
singlet  lies  lower  as  one  can  see  from  Table  3-3.  Of  those  states  in  which  one  electron 
is  in  each  band  (F&^  type),  it  is  clear  that  the  lowest  state  of  A2  symmetry  is  a triplet 
lower  than  any  other  state  of  the  F ^ type.  Of  the  wave  functions  in  which  both  electrons 
are  in  the  same  band  (F  type),  the  lowest  state  is  a singlet  of  B.  symmetry.  As  to 
the  competition  between  these  two,  the  lowest  one  will  be  determined  by  the  actual 
values  of  the  parameters  involved.  Either  one  of  these  two  could  lie  lower.  Thus  we 
see  that  in  the  case  of  non-degenerate  bands  the  singlet  lies  lower  whereas  in  the  case 
of  degenerate  bands  it  is  possible  that  the  triplet  may  lie  lower. 
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(A  TWO- ELECTRON  EXAMPLE  OF  FERROMAGNETISM) 


Energies  for  two-electron  states  arising  from 
degenerate  bands  in  a square  lattice 


1 

II 

in 

IV  v 

singlet 

1 

I 

E4-(VoO)EJ-U(<r/.  <f/)E2.8(Vo.  gf)  E 

♦64(  g*  - <f,V*0 

E ■ V ♦J... 
o 

• £ *(l  -cV 

E. ! 

(Vo«  J)(.  ♦ci) 

E.t  lji./gl*  tf-gf — T-  - 

' 1 1 (VQe  J)(l  ♦ c2) 

singlet 

3 

E 

E ■ 0 

E » 0 

singlet 

Bj 

1 

E4-(Vo-J)EJ-I*(<;'12.  <;2Z)E2  * *(Vo  - J){«f  ^ ♦ g/) E 

♦ <J2V-  o 

E ■ V -J... 
o 

»<?/(!  -cV 

E- ' j-  ... 

(V  - J)(l  *cZ) 
o 

e •tiyijg**  g,*- g*  l6c*  ... 

‘ 1 ' (Vo.J)(|.c*) 

singlet 

B 

E * 0 

E-0 

triplet 

E 

Fss 

Eatiyitf, 

E.ti,/^ 

E*t  l/l  g^ 

E*4-  lVigi 

singlet 

B 

F.b 

E-0 

E ■ 0 

singlet 

B 

F.b 

E "0 

E-0 

triplet 

F.b 

E,-E<(K-J)-8(«f1«  <f2)EM(<f,*  (f^lK -J).  0 

E a K - J*. 

E. gj.-f-ig^g/,... 

singlet 

3 

F.b 

E = 0 

E « 0 

triplet 

3 

F.b 

E *0 

E * 0 

singlet 

F.b 

EJ-E*(K«  J)-8 (<?,  ♦ <fj)E«4(  J)a0 

E * K ♦ J ♦ . . . 

singlet 

F.b 

E-t  Silgx*gt) 

e * t Ji(  gx  ♦ gj) 

triplet 

E 

F.b 

E.tyi(<f,*  <?2> 

E-t  Sligx<  gz) 
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4.  A SYSTEM  OF  TWO  p ELECTRONS 


In  connection  with  the  work  on  the  two-electron  example  of  ferromagnetism 
developed  by  other  members  of  the  Group,  ^ there  occurred  a question  regarding  the 
energy  levels  of  a system  of  two  Ni  atoms  as  a function  of  the  internuclear  distance; 
the  reason  is  that  the  properties  of  the  pair  of  atoms  might  be  supposed  to  give  infor- 
mation about  the  ferromagnetic  substance  in  bulk. 

In  order  to  obtain  qualitati  /e  information  about  this  system,  particularly 
about  the  multiplicity  of  the  lowest  lying  state,  the  problem  was  simplified  several 
steps  as  follows.  A Ni  atom  which  has  one  hole  is  3d  orbits  is  first  reduced  to  an  atom 
with  a hole  in  p orbits,  because  both  d and  p levels  give  rise  to  degenerate  bands. 

A hole  in  p orbits  behaves  qualitatively  similar  to  an  electron  in  p orbits.  Further- 
more, the  Gaussian  type  of  wave  functions  may  be  enough  to  give  the  qualitative  infor- 
mation. 

After  these  considerations,  now  the  energy  levels  of  a system  of  two  atoms 
each  having  one  p electron  is  being  calculated  using  Gaussian  wave  functions,  as  a 
function  of  internuclear  distance,  and  taking  into  account  the  configuration  interaction. 
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5.  SPIN  OPERATOR  METHOD  --  APPLICATION  TO  CONFIGURATION  INTERACTION 


A detailed  account  of  an  operator  method  of  constructing  eigenfunctions  of  S2 
was  presented  in  the  last  Progress  Report.  The  application  of  this  technique  to  con- 
figuration interaction  has  proved  rather  successful  and  some  of  the  details  of  its  use 
will  be  given  here.  The  following  remarks  will  be  largely  confined  to  the  treatment  of 
the  Mn~  -Q^ — Mn++  model  of  antiferromagnetism  discussed  elsewhere  in  the  present 
Progress  Report. 

It  was  shown  in  the  April  15  Quarterly  Progress  Report  that  there  exists  an 
operator  O when  operating  on  a 2N  electron  determinant  whose  value  is  zero, 
produces  a linear  combination  of  2N  electron  determinants  such  that  the  resulting  state 
is  a singlet.  The  only  condition  imposed  was  that  all  the  one -electron  spatial  orbitals 
be  distinct,  but  not  necessarily  orthonormal.  It  was  also  pointed  out  in  that  Report 
that  this  operator  method  could  be  used  to  set  up  the  matrix  components  of  the  Hamil- 
tonian between  singlets  corresponding  to  different  configurations.  Suppose  Pj  is  a 
product  function  of  space -spin  orbitals 

Pl  = Ujd)  a(  1 ) u2(2)  a(2)  ...  uN(N)  a(N)  uN  + j(N  + 1)  p(N  + l)  ...  u2N(2N)  p(2N)  (5.  1) 

of  total  Mg  = 0.  Then  one  of  the  singlet  states  that  can  be  formed  from  this  collection 
of  spatial  orbitals  is 


4-  = AOpj  (5.2) 

where  O is  the  spin  operator  and  A is  the  antisymmetrizing  operator.  The  matrix 
component  of  the  energy  between  two  configurations  for  this  particular  type  of  singlet 
is 


<«!»,  H +2>  = <AOpj  H AOp2>  (5.3) 

We  now  make  use  of  the  following  properties  of  the  spin  operator  O:  (1)0  commutes 
with  a spin-free  Hamiltonian.  (2)  O commutes  with  the  antisymmetrizing  operator  A. 
(3)  The  square  of  O is  equal  to  a constant  times  O.  These  properties  allow  one  to 
write  Eq.  (5.  3)  as 

<4'1IHI4<2>  * (N  + l)1/2  <A(2N)  0(2N)  Pj  IHI p2  > (5.4) 

where  pj  and  p2  are  spin  product  functions  corresponding  to  different  configurations 
and  A^N)  and  signify  that  these  operators  are  functions  of  the  coordinates  of  all 

of  the  electrons. 

In  the  ground  configuration  of  the  Mn^  - — Mn++  system  the  orbitals  con- 
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sidered  are  the  five  Mn  3d  orbitals  each  singly  occupied  and  the  oxygen  2p  orbitals 
which  arc  assumed  to  hold  six  electrons.  Thus  we  have  a sixteen-electron  system  but 
only  thirteen  distinct  orbitals.  The  object  of  these  remarks  is  to  generalize  the  spin 
operator  O in  order  to  handle  systems  composed  of  spatial  orbitals  all  of  which  are 
not  distinct.  By  doing  this  we  will  be  able  to  cast  matrix  elements  in  the  convenient 
form  given  in  Eq.  (5.  4). 

In  the  investigation  of  the  Mn^^— — Mn++  system  the  singlet  considered  for 
the  ground  configuration  was  formed  by  combining  two  S = 5/2  states  made  up  from  the 
Mn  3d  orbitals  to  give  a resultant  ten-electron  singlet.  The  remaining  six  electrons 
were  assigned  to  the  three  2p  oxygen  orbitals.  In  terms  of  the  spin  operator  this  six- 
teen-electron singlet  is  written 

$ = A(l6)  [uj(1)  o(  1)  Uj(2)  0(2)  . . . u3(5)  q(5)  u3(6)  0(6)  • O(l0)  {u4(7)  a(?)  u&(8)  a(8) 

ii  (5.  5) 

u6(9)  a(9)  ...  ug(ll)  a(  1 1)  u9(12)  0(  12)  u1Q(l3)  0(13)  ...  un(l6)  0(l6)JJ 

where  u j,  u^,  and  u3  represent  the  2p  functions  and  u^  through  Ug  are  the  3d  functions 

on  one  Mn++  ion  and  uQ  through  u,  , are  the  3d  orbitals  on  the  other  Mn++  ion.  An  es- 

9 1 J (10) 

sential  point  appears  here;  that  is  the  ten-electron  spin  operator  O'  ' does  not  com- 
mute with  the  sixteen-election  antisymmetrizing  operator  A^l6\  Therefore,  unless 
this  singlet  can  be  set  up  using  the  full  sixteen-electron  spin  operator  the  ma- 

trix components  cannot  be  ha  idled  as  in  (5.  4). 

Consider  a 2N  electron  system  with  2(N  - n)  distinct  orbitals  which  will  be 
kept  singly  occupied  and  n distinct  orbitals  which  will  be  doubly  occupied.  The  or- 
bitals may  be  catalogued  in  the  following  manner 


Group  1^ 


“N, 


n+  1, 


Group  11^ 


u 


B 


Group  Ig 


Group  Ilg 


(5.6) 


If  all  of  the  functions  were  distinct,  the  2N  electron  singlet  under  discussion  would  be 
written 


* = aUN)q(2N)  |u  (l)  a(  1)  ...  uN^(N)  a(N)  u^N*  1)  0(N+  1)  . ..  u^UN)  0(2N)}(5.  7) 
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Suppose  now  we  set  Uj^  equal  to  uj^  for  j going  from  1 to  n in  the  singlet 
given  in  (5.  7).  Such  a change  cannot  affect  the  multiplicity  of  the  total  wave  function. 
However,  many  of  the  determinants  which  comprise  (5.  7)  will  vanish  as  a result  of 
assigning  the  same  spin  to  orbitals  whose  spatial  parts  are  now  assumed  identical.  Of 
the  non-vanishing  determinants,  many  will  differ  from  others  merely  by  interchanges 
of  rows  and  hence  these  determinants  can  be  combined.  The  net  result  of  this  change 
of  spatial  orbitals  is  that  the  normalization  of  the  resultant  singlet  is  destroyed.  O 
must  be  modified  so  that  the  normalization  can  be  preserved  under  such  a change. 

It  can  be  shown  that  the  necessary  modification  of  is  to  multiply  it  by  a 

constant  factor  which  depends  only  on  the  total  number  of  electrons,  2N,  and  the  num- 
be  of  identical  orbitals,  n.  This  factor  is 


£ n*  (N  - k)! 
k = o<n  - *);  N! 


1/2 


(5.8) 


Once  the  matrix  elements  have  been  reduced  to  the  form  of  Eq.  (5.  4),  O^2*^ 
can  be  replaced  by  an  equivalent  operator  for  the  purpose  of  taking  the  matrix  com- 
ponent. The  equivalent  operator  is  obtained  by  omitting  entirely  the  one -electron  op- 
erators Sj  , S/  where  i and  j refer  to  electrons  occupying  an  orbital  which  is  doubly 
filled.  More  precisely,  this  operator  is  O2^  ” and  it  operates  only  on  the  spin  co- 
ordinates of  the  electrons  in  groups  and  Ilg.  This  is  a valuable  result  as  it  ma- 
terially reduces  the  amount  of  work  necessary  to  find  the  matrix  components.  There- 
fore, in  the  case  of  Mn — O — Mn  the  matrix  component  of  the  energy  between  the  singlets 
corresponding  to  the  i**1  and  j**1  configurations  can  be  expressed  as 


<*llHl*j>-  (N  - n + 1)1/2<A(16)  O(10)  PiIH|Pj> 


(5.  9) 


where  N is  8 and  n is  equal  to  the  number  of  doubly  filled  orbitals  which  is  three. 


G.  W.  Pratt,  Jr. 
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6.  ANTIFERROMAGNETISM 


At  the  present  time  attempts  to  give  a quantum  mechanical  description  of  anti- 
ferromagnetism have  been  largely  restricted  to  the  treatment  of  simple  models  which 
supposedly  would  show  a preference  for  spins  to  be  oriented  antiferromagnetically.  A 
configuration  interaction  treatment  of  the  Mn*-*— — Mn**  system  is  reported  here. 

The  Mn — Mn  distance  as  determined  by  x-ray  scattering  is  4.  426  A0^  and  it  seems 
unlikely  that  the  Mn  3d  orbitals  located  about  centers  this  distance  apart  would  overlap 
appreciably.  This  overlap  was  neglected  in  this  investigation.  By  considering  excited 
configurations  in  which  there  is  an  unpaired  spin  on  the  oxygen,  a spin  dependent  coupl- 
ing will  appear  which  causes  the  non-overlapping  charge  distributions  associated  with 
the  Mn  ions  to  interact  through  the  oxygen.  This  is  commonly  referred  to  as  super- 
exchange interaction.  ^ The  object  of  this  work  was  to  examine  the  nature  of  this 
coupling  to  determine  whether  it  would  tend  to  orient  the  spins  of  the  Mn  electrons 
ferromagnetically  or  non-magnetically. 

The  ground  configuration  was  taken  to  be  two  Mn  * ions  and  an  O-  ion.  The 
Mn**  ion  has  five  3d  electrons  outside  closed  shells  and  the  O-  ion  was  taken  to  have  a 
rare  gas  configuration.  The  excited  configurations  considered  were  those  obtained  by 
allowing  an  electron  in  an  orbital  associated  with  the  0~  ion  to  occupy  one  of  the  Mn  3d 
orbitals.  The  only  symmetry  considered  in  the  work  was  of  the 
state  has  this  symmetry. 

The  one-eleciron  functions  used  were  assumed  to  be  localized  and  orthonormal 
although  the  results  obtained  are  independent  of  their  exact  nature.  Table  6-1  shows 
the  assignment  of  electrons  for  the  various  configurations.  The  Mn  3d  orbitals  are  de- 
noted by:  mj  = 2 (6+),  m ^ = 1 (*+),  m j = 0 (v),  nij  = - 1 (»_),  and  mj  = - 2 ( 6 _ ) . 

The  oxygen  orbitals  are  m } = 1 (P+),  nij  = 0 (PQ),  “j  E -1  (P_). 

The  states  of  maximum  multiplicity,  S = 5,  for  each  configuration  and  the 
corresponding  secular  equation  was  set  up.  The  energy  of  the  lowest  S = 5 state  was 
compared  with  the  energy  of  the  lowest  singlet  found  by  constructing  one  particular 
singlet  for  each  configuration  and  examining  the  resulting  secular  equation.  Let  us  con- 
sider the  choice  of  this  singlet  state. 

For  each  configuration  there  are  ten  singly  occupied  orbitals  and,  therefore, 

42  orthogonal  singlets.  From  neutron  diffraction  studies  of  MnO,  it  is  thought  that  the 
Mn**  ions  can  be  divided  into  two  sublattices  each  sublatHce  being  ferromagnetically 
coupled:^  the  singlet  state  most  closely  resembling  this  situation  is  that  formed  by 
combining  the  two  sublattices  which  are  separately  in  a state  of  maximum  spin  to  give  a 
resultant  singlet  for  the  entire  system.  The  analogue  of  this  picture  in  the  model  being 
considered  here  is  for  the  two  Mn  ions  to  be  each  in  a state  of  maximum  multiplicity 
and  for  these  two  S = 5/2  systems  to  combine  to  give  a singlet  ftr  the  ground  configura- 
tion. The  singlet  for  each  excited  configuration  was  formed  by  the  combination  of  two 
S = 5/2  systems.  As  it  is  only  possible  to  form  one  independent  singlet  in  this  manner. 


T*  type  as  the  ground 

S 
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Table  6-1 


for  each  configuration  one  S = 5 and  one  S = 0 state  was  considered. 

The  states  which  are  even  with  respect  to  a reflection  in  a plane  perpendicular 
to  the  symmetry  axis  and  passing  through  the  oxygen  are: 


*1*  ■ +2  ♦ +3 
#2g  - *4  * +5 
*3*  ' +6  * +7 

•«* 

*5*  = +10  * +11 


From  these  six  states  the  states  which  have  the  + symmetry  with  respect  to  reflection 
in  a plane  containing  the  symmetry  axis  of  the  model  can  be  formed.  These  four  states 
are: 


*1  = #o8  = +1 

*2  #lg  + #4g  = *2  + *3  + +8  + *9 
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*3  = #28  + #3g  = *4  + +5  + *6  + *7 
♦4  = #5g  = *10  + *11 


(6.  2)  con'd. 


The  S = 5 and  S = 0 states  corresponding  to  these  four  £+  states  were  set  up 
and  the  corresponding  secular  equations  written  down.  The  S = 5 states  were  handled 
by  conventional  methods  in  that  the  states  ^(S  = 5)  are  just  a single  determinant. 

The  singlet  states  4^(S  = 0)  are,  however,  much  more  complicated  as  these  singlets 
are  each  expressed  as  a linear  combination  of  252  determinants.  It  would  be  an  enor- 
mous task  to  compute  the  matrix  components  of  the  energy  for  these  states  by  the  cus- 
tomary procedures.  These  matrix  components  were  found  by  the  spin  operator  method 
reported  in  the  preceding  Quarterly  Progress  Report.  The  essence  of  this  method  is 
that  one  can  construct  a state  of  definite  multiplicity  by  operating  on  a single  determi- 
nant of  one -electron  space -spin  orbitals  with  an  operator  which  is  essentially  a poly- 
nomial the  elements  of  which  are  step-up  and  step-down  spin  operators  taken  to  vari- 
ous powers.  This  operator  has  the  properties  that  it  is  hermetian,  it  commutes  with 
a spin  free  Hamiltonian,  it  commutes  with  the  S*  operator,  it  commutes  with  the  anti- 
symmetrizing  operator,  and  its  square  is  equal  to  a constant  times  itself.  These  prop- 
erties allow  one  to  use  the  operator  to  take  matrix  components  of  the  Hamiltonain  be- 
tween singlets  corresponding  to  all  configurations.  ^ Some  of  the  details  involved  are 
given  in  another  section  of  the  current  Progress  Report. 

In  terms  of  the  spin  operator  O and  the  antisymmetrizing  operator  A,  we 
can  write  the  four  1 Y + states  as 


♦ i(S  « 0) 

♦2(S  = 0) 

♦3(S  = 0) 

♦4(s  = 0)  = ACXpl0  + Pu) 


A Op  j 

AO(p2  + p3  + Pg  + p9) 
ACXp4  + p5  + P6  + P7) 


(6.  3) 


where  the  p^  are  spin  product  functions  referring  to  the  various  configurations. 


Pj  = 6+(l)  a(  1)  6 (2)  a(2)  x+(3)  a(3)  w_(4)  a(4)  <r(5)  a(5)  Pq(6)  a(6)  Pq(7)  0(7)  P+(8)  0(8)  P+(9)  0(9) 

P_(10)  a(  10)  P_(  1 1)  0(  1 1)  w*(  12)  0(  12)  w_*(l3)  0(  13)  w+'(l4)  0(  14)  6_'(  15)  0(  15)  6+  '(16)  0(l6) 

(6.4) 
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p2  * «+( 1)  a(  1)  . . . d 5)  a(5)  PQ( 6)  q(6)  Pq(7)  0(7)  P+(8)  a(8)  P+(9)  0(9)  P_(  10)  a(  10)  5 ( 1 1)  0(  1 1) 
<r'(l2)  0(12)  . . . 6+'(l6)P(l6)  (6.5) 

P3  * *+(  1)  a(  1)  . . . «r(5)  a(5)  Pq(6)  a(6)  Pq(7)  0(7)  P+(8)  a(8)  P+(9)  0(9)  6_ '( 1 0)  a(  10)  P_(  1 1 ) 0(  1 1) 
.'(12)0(12)  ...  6+ '(16)  0(16)  (6.6) 

P4  = 6+(  1)  a(  1)  . . . <r<5)  a(5)  Po(6)  a(6)  Pq(7)  0(7)  P+(8)  a(8)  P+(9)  0(9)  P_(10)  a(l0)  w_(  1 1)  0(l  1) 
.'(12)0(12)  ...  6+  '(16)0(16)  (6.7) 

P5  = 6+(l)a(l)  ...  <r(5)  a(5)  Pq(6)  a(6)  PQ(7)  0(7)  P+(8)  a(8)  P+(9)  0(9)tr_'(lO)a(lO)P_(ll)0(ll) 
.'(12)0(12)  ...  6+ '(16)0(16)  (6.8) 

p6  = 6+(l)a(l)  ...  ^5)Q(5)Po(6)o(6)Po(7)0(7)P+(8)a(8)«+(9)0(9)P_(lO)a(lO)P_(ll)  0(ll) 
.'(12)0(12)  ...  6+ '(16)0(16)  (6.9) 

p?  = 6+(  1)  a(l)  . . . «r(5)  a(5)  Pq(6)  q(6)  Pq(7)  0(7)  w+'(8)  a(8)  P+(9)  0(9)  P_(  10)a(  10)  P_(  11)0(11) 
.'(12)0(12)...  6+'(l6)0(l6)  (6.10) 

pg  = 4+(l)a(l)  ...  <r(5)Q(5)Po(6)a(6)Po(7)0(7)p+(8)a(8)6+(9)0(9)P.(lO)Q(lO)P_(ll)0(ll) 
.'(12)0(12)  ...  6+ '(16)0(16)  (6.  11) 

P9  = 6+(  1)  a(  1)  . . . <r(5)  a(5)  Pq(6)  a(6)  Pq(7)  0(7)  6+'(8)  a(8)  P+(9)  0(9)  P_(  10)  a(  10)  P (11)  0(11) 
.'(12)0(12)...  6+  '(16)0(16)  (6.12) 

p10  = 6+(l)a(l)  ...  ^5)a(5)Po(6)a(6)«r(7)0(7)P+(8)Q(8)P+(9)0(9)P_(lO)a(lO)P_(ll)0(ll) 
.'(12)0(12)  ...  6+'(l6)  0(16)  (6.13) 

Pn  « 6+(l)a(l)  ...  <K5)a(5)  a*(6)  a(6)  Pq(7)  0(7)  P+(8)  a(8)  P+(9)  0(9)  P_(10)  o(l0)  P_(l  1)  0(1 1) 

.'(12)0(12)  ...  6+'(l6) 0(16)  (6.  14) 
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The  matrix  elements  of  the  energy  for  the  S = 0 case  are  given  as: 


H11 

= ^AOpjIHIpj) 

(6.15) 

H 12 

= y(AOp1IH|p2  + P3  + p8  + p9) 

(6. 16) 

H13 

= y(AOp1lH|p4  + p5  + p6  + p?) 

(6.  17) 

H14 

= Y(AOpjlH|p10  + pn) 

(6.  18) 

H22 

= y(AO  jp2  + p3  + pg  + P9J  IH|  jp2  ♦ p3  ♦ pg  ♦ p9}  ) 

(6.19) 

H23 

* y(AO  |p2  + P3  + Pg  + P9}  IHI  {p4  + P5  + p6  + p7}  ) 

(6.20) 

H24 

= y(AO  |p2  + p,  ♦ p8  + p9}  IHI  [p10  ♦ pu}) 

(6.21) 

H33 

= y(AO  jp4  + P5  ♦ p6  + p7}  IH|  |p4  + p5  + p6  + p?}  ) 

(6.22) 

H34 

= y(AO  |p4  + p5  + p6  + p7}  IHI  {p10  ♦ Pj  i } ) 

(6.  23) 

H44 

= Y(AO  {p10  + pu}  IHI  {p10  + Pj  i}) 

(6.24) 

Y = Vb 

As  an  illustration  of  how  the  spin  operator  method  works  the  element 
y(AOp2/H/p^)  will  be  worked  out  in  detail.  All  of  the  one-electron  orbitals  are  ortho- 
normal  and  it  has  been  assumed  that  there  is  no  overlap  between  the  orbitals  localized 
about  the  two  Mn  ions.  The  spin  operator  used  can  be  expressed  as 

O = i (l  ♦ 

Y t 

N N N 

ci  ^ 's;X  > * c2  l <s;A k.  sL sL>  * ■ ■ ■} 

JAkB  A B JA>kAmB>nB  A A B *B 

(6.25) 

where 

„ _ ( nM  (N  - M)i 

' N!  M! 

Since  O is  to  operate  on  a spin  product  function,  we  may  speak  of  the  electrons  as 
being  distinguishable.  Electrons  1,  3,  4,  5,  and  10  comprise  set  A.  Electrons  12, 
13,  14,  15,  and  16  make  up  set  B.  N is  the  number  of  electrons  in  set  A or  B and 
M is  the  number  of  {J  spins  assigned  to  the  electrons  in  set  A or  the  number  of  a 
spins  assigned  to  set  B.  Therefore,  the  term 
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C 


1 


N 

z 

JAkB 


in  O means  that  runs  over  the  electrons  in  set  A and  kg  is  to  run  over  the  elec- 
trons in  set  B.  (AOp^/H/p^)  written  out  is 


AO  {6+(  1)  a(  1)  6 (2)  a(2)  w+(3)  a(3)  w_(4)  a(4)  <r(5)  a(5)  Pq(6)  a(6)  Pq(7)  p(7)  P+(8)  a(8) 
P+(9)p(9)P_(l0)a(l0)6_(ll)p(ll)<r'(l2)p(l2)ir_'(l3)p(l3)*+'(l4)p(14)6_'(15) 
P(l5)6+'(l6)p(l6)}  I HI  {6+(l)a(l)6_(2)a(2)w+(3)a(3)*_(4)Q(4)a(5)a(5)Po(6)  (6.26) 

a(6)  Pq(7)  p(7)  P+(8)  a(8)  P+(9)  P(9)  »_'( 10)a(  10)  P ( 1 1 ) p(  1 1)  <r'(  12)  p(  12)  ir_'(  1 3) 
P(l3)w+'(14)P(14)6_'(15)  P(15)6+'(16)P(16)} 


where  the  orbitals  which  do  not  match  are  underlined. 

The  1 from  0 combined  with  the  1 from  A gives  for  a contribution 

(P_(  10)  6 (11)  | ir_ '( 10)  P ( 1 1)) 


Due  to  the  spatial  orthogonality  of  the  orbitals  and  due  to  spin  orthogonality,  the  only 
contributing  permutation  from  A that  can  be  combined  with  the  1 from  0 is  the 
identity  permutation. 

The  term 


C 


1 


N 

I s- 

jAkB  A 


combined  with  the  1 from  A makes  no  contribution  due  to  spin  orthogonality.  The 
only  term  to  give  a non-zero  result  is  Cj  S~q  Sjj  which  must  be  combined  with  P^q 
P^O  j ^ from  A,  the  combination  yielding 

- C j(6_(  1 1 ) P_(  13)  f P_(  11)  w_'(  1 3)) 

There  will  be  no  other  terms  entering  into  this  matrix  element.  Therefore, 

(AOp2)IHIp5)  = (1  - Cj)(6_(l)  P_(2)IP_(  1)  w_'(2)  (6.  27) 
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and  Cj  is  - 1/5.  By  proceeding  in  exactly  the  same  manner,  all  of  the  matrix  elements 
of  H for  the  S = 0 case  can  be  set  up. 

Let  the  Hamiltonian  matrix  for  the  singlet  states  be  denoted  by  Hq  and  the 


Hamiltonian  for  the  S = 5 states  be  H^.  Then  we  can  set 

HQ  = H5  + 


The  matrix  D has  the  following  form 


(6.  28) 


o 

o 

o 

o 

o 

°22 

°23 

o 

o 

°23 

°33 

o 

o 

O 

O 

D. 

(6.29) 


Now  suppose  4»0  is  the  eigenvector  of  HQ  corresponding  to  the  lowest  eigen- 
value of  HQ  which  we  shall  call  Eg.  Then 

(4»0IHl4<0)  = (4.0IH5I*0)  + (+0lDl4»0)  (6.30) 


E0  = * 0,  5 * D 


If  4»5  is  the  eigenvector  of  H^  corresponding  to  the  lowest  eigenvalue  E^,  then  by  the 
variational  principle 


(*5IH5I*5)  $ (*0IH5I4»0) 


(6.31) 


E5  * *0.  5 


Therefore, 


Eq  > E5  + D 


(6.32) 
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W-  f j«  * "t*  '**  **•'  ?>-•  V “-«  8 respectively,  and  for  k = 12, 

■|jg  fl^jl.6.  ■■  •*'• '»  »Wl  6,‘  respectively.  %3  is  obtained  by  replacing 

-£"-  _i/y  ajjj  j^»  by  »^«.  everyv^Itt’i’U-in  O^,  except  when  4>+  and  &+ 1 occur  as  the  result 


^ *U(*)):  + *6 J^UjlPjd)  tr_'(Z))  * (*_(>>  P+(2)iP+(0  Tt_'(2)) 

b|(6  ’(DP  (2)|  P (1)  rr  (2))  + <6  ’( l)  P (2)|P  ( D * (2))  +•  (6  '( l)  P , <2)|P  A 1 ) tr  (2)) 

■>  |(6T(r)  Pt(2)|Pf(l)Trr,(2))  + (6r(D  P0U)|P(j(l)7rT'(2))  i (6+(l)  P_(2)|  P_(  1 ) tt+ '(2)) 

+ |(6+.'a)P+(2)|P+U)ir+(2))  c (6/(  DP0(2)|P_U)trtl2f)  + (6+l.(l)  P_-(2)|P-  ( D tr+(2)) 

(6.  34) 


+ 


|kZ2(Vu)p0w|p0(.)+K'(k)  * |ki  (♦,(') 

-yt<r(  l ) PoU)|P0<  1 ) « '(2))  ’ + 2(a(  1 ) P+(2)|P+i  lj  o 


Po(k)|PolD<)*k(k)) 

'(2))  c 2(tf(l)P_(2)|P  U)<rK2)) 


(6.  35) 

Tfitr principal  parts  of  and  are  exchange  integrals  between  the  oxygen 
orbitals  and  the  Mn  orbitals;  these  integrals  are  all  positive.  The  other  integrals  en- 
tering D.? , and  are  three-center  integrals  which  cannot  be  of  major  import  as  they 
involve  the  orbitals  of  both  Mn  ions  Although  these  three-center  integrals  are  by.  no 
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A 16  A 5 

D22  = | I («k'(l)P.(k).'P.(l)«k,(k))  < | I (♦k(l)P»|P.(l)«k(k)) 

k “ 12  k - 1 

16  A 5 

+ | I (♦k’(l)P+(k)IP+(l)«k'(k))  + | Z (♦k(l)P+(k)IP+(lHk(k)) 
k * 12  k - 1 

♦ H(P_(1)6_U)|6_'(1)P_U))  + 2(6_(1)P0(2)IPo(1)&_'(2))  + 2(6_(l)  P+(2)|P+(l)  6_'U)) 

+ -^(P+(l)6+(2)|6+’(l)P+(2))  + 2(6+(l)Po(2)|Po(l)6+'(2))  + 2(6+(  1)  P.(2)|P_(  1)  &+ '(2)) 

(6.  33) 

We  mean  by  <j>k  fork  = 1.  2,  3,  4,  5;  6+,  6 , w+,  w_,  a respectively,  and  for  k = 12, 
13,  14,  15,  1 6;  cr',  w_\  w+\  6 and  respectively.  D33  is  obtained  by  replacing 
6+  by  w+  and  6+'  by  w+'  everywhere  in  D^2  except  when  6+  and  6+'  occur  as  the  result 
of  a summation. 

°23  = |(6.(1)P.(2)|P_(1)*_'(2))  + (6_(l)Po(2)|Po(l)ir_’(2))  + (6_(l)  P+(2)|P+(  1)  w_'(2)) 

+ |(6_'(l)P_(2)|P_(l)w_(2))  + (6_,(l)P0(2)|Po(l)n  (2))  + (&  ’(l)  P+(2)|P+(l)  ir_(2)) 

+ |(6+(l)P+(2)|P+(l)w+’(2))  +(6+(l)P0(2)|Po(l)w+-(2))  + (&+(  1 ) P_(2)l P_(  1 ) w+  '(2)) 

+ |(6+'(l)P+(2)|P+(l)w+(2))  +(6+’(1)P0(2)|P0(1)it+(2))  + (6+'(l)  P_(2)|P_(l)  w+(2)) 

(6.  34) 

16  5 

D«  = 5k£2(v(*>po(k»po“>v(k>  * 

+ -^Ml)P0(2)|Po(l)tr'(2))  + 2(a(l)P+(2)jP+(l)a’(2))  + 2(«r(l)  P_(2)|P_(l)  tr'(2)) 

(6.  35) 

The  principal  parts  of  and  D33  are  exchange  integrals  between  the  oxygen 
orbitals  and  the  Mn  orbitals;  these  integrals  are  all  positive.  The  other  integrals  en- 
tering D and  Dj3  are  three-center  integrals  which  cannot  be  of  major  import  as  they 
involve  the  orbitals  of  both  Mn  ions.  Although  these  three-center  integrals  are  by  no 
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means  zero,  they  ? re  overshadowed  by  the  exchange  Integrals.  Therefore  and 
D33  as  well  as  D44  are  positive.  Since  D44  is  positive,  if  D is  to  have  a negative 
eigenvalue  it  must  come  from  the  2 x 2 secular  equation  involving  D,,,  D,,,  and  D,v 
The  secular  polynomial  coming  from  this  part  of  D is 

x2  * {D22  + °33>  X + D22D33  ‘ <D23>2  = 0 (*•  W 

We  can  only  have  a negative  root  if  D22D33  < ^°23^  since  the  sum  of  the  roots  is 
positive.  By  means  of  the  same  arguments  which  established  the  positive  character 
of  D22  and  D33  we  can  see  that  D23  is  less  than  either  D22  or  D3J.  Therefore,  the  D 
matrix  has  no  negative  roots  and  consequently  the  energy  of  the  lowest  lying  singlet 
state  as  found  by  this  configuration  interaction  cannot  lie  below  the  lowest  S = 5 state. 

Up  to  this  point  configurations  in  which  an  electron  on  the  oxygen  ion  goes  into 
a 4s  orbital  on  one  of  the  Mn  ions  have  not  been  included.  The  4s  orbitals  are  more 
extended  than  the  3d's  and  now  we  would  have  to  admit  exchange  integrals  between  a 
4s  orbital  localized  about  one  Mn  ion  and  the  3d  orbitals  localized  about  the  other  Mn. 

By  examining  the  character  of  D22  and  D33  it  appears  that  the  inclusion  of  these  con- 
figurations could  not  lead  to  a lower  lying  singlet  energy  than  that  of  the  S = 5 states. 

There  is  one  quite  definite  conclusion  which  can  be  drawn  from  this  investiga- 
tion. Within  the  framework  of  the  approximations  made  here  it  seems  clear  that  when 
an  electron  is  removed  from  the  O”  ion,  leaving  an  unpaired  spin  on  the  oxygen,  that 
this  unpaired  spin  couples  the  non-overlapping  charge  densities  on  the  Mn  ions  ferro- 
magnetically.  In  view  of  this  result,  it  appears  likely  that  some  other  appreaen  would 
be  more  profitable  as  an  explanation  of  antiferromagnetism.  Of  course,  it  may  well 
be  that  the  singlet  state  considered  here  in  no  way  resembles  the  antiferromagnetic 
state. 

Before  abandoning  this  simplified  model  there  is  one  feature  of  the  configura- 
tion interaction  treatment  that  bears  examination.  That  is  the  fact  that  this  approach 
has  placed  all  of  the  emphasis  on  the  probability  amplitudes  of  the  various  configura- 
tions, the  energy  being  stationary  with  respect  to  the  variation  of  these  parameters. 

It  is  quite  possible  that  the  one -electron  wave  functions  are  really  quite  different  for 
the  singlet  and  S = 5 states.  Furthermore,  the  use  of  orthogonal  orbitals  may  force 
one  to  carry  out  a rather  extensive  configuration  interaction  to  accurately  describe  the 
correct  state  of  affairs.  An  investigation  is  being  carried  out  at  present  using  a mix- 
ture of  orthogonal  and  non -orthogonal  orbitals  but  considering  only  the  ground  configura- 
tion. The  energy  will  be  stationary  with  respect  to  variations  of  parameters  in  the  one- 
electron  functions  and  these  functions  will  be  determined  separately  for  the  states  of 
different  multiplicity.  Preliminary  results  obtained  for  this  approach  have  been  en- 
couraging. 
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7.  AN  EXTENSION  OF  HUND  'S  RULE 


In  the  course  of  configuration  interaction  calculations  with  orthogonal  orbitals 
it  is  sometimes  interesting  to  know  what  the  multiplicity  is  of  the  lowest  state  of  some 
given  class  of  states.  In  general  it  is  impossible  to  answer  this  question  without  cal- 
culations. Under  one  set  of  circumstances  it  is  possible  to  answer  this  question  with- 
out any  calculation.  We  know  of  Hund's  rule  of  atomic  spectra.  Part  of  this  rule  states 
that  of  all  the  levels  arising  from  a given  spacial  configuration  the  state  of  highest  mul- 
tiplicity lies  lowest.  A rule  similar  to  this  can  be  proved  in  general. 

We  imagine  n distinct  orthogonal  orbitals  u^,  u^,  ...  un.  We  put  one  elec- 
tron in  each  of  these  orbitals  with  either  spin  up  or  spin  down.  This  gives  rise  to  2n 
states.  These  states  can  be  formed  into  states  of  definite  multiplicity  ranging  from 
singlets  (for  n even)  to  states  of  multiplicity  n + 1.  (The  presence  of  other  orbitals 
which  are  doubly  filled  does  not  influence  the  validity  of  the  following  arguments  and 
can  be  dropped  from  considei  ition.  ) We  shall  now  show  that  there  is  no  state  that  lies 
lower  than  the  state  of  highest  multiplicity.  We  lose  no  generality  by  restricting  our- 
selves to  all  those  states  with  the  same  z component  of  total  spin.  Of  these  only  one 
corresponds  to  the  state  of  highest  multiplicity  whereas  for  any  other  multiplicity 
there  is  in  general  more  than  one  state  of  this  multiplicity.  Since  states  of  different 
multiplicities  cannot  interact  we  are  confronted  with  the  problem  of  diagonalizing  for 
each  multiplicity  the  matrix  of  the  Hamiltonian 

H = £ - V*  ♦ f(i)  + I (7.  1) 

i > j J 

Here  f(i)  is  the  one -electron  potential  and  g. . is  the  interelectronic  Coulomb  interac- 

( 1 ) 

tion.  From  the  Dirac  vector  model  we  know  that  for  a given  multiplicity  the  matrix 
of  the  Hamiltonian  can  be  written  down  very  simply, 

*\*  = C°nSt-  + £ • Jij(Pij)V<r  (7-2) 

ij 

Here  represents  the  Hamiltonian  matrix  and  (Py)^  are  unitary  matrices  which 

represent  the  permutation  of  the  spatial  coordinates  of  electron  i with  those  of  elec- 
tron j.  Jy  is  the  customary  exchange  interaction  Jy  = Juf(l)  u^*(l)  g^  u^U)  u^(2) 
dTi  d-r^  which  is  greater  than  zero  since  it  represents  the  interaction  of  a charge  den- 
sity with  itself.  We  first  consider  the  state  of  highest  multiplicity.  This  has  a Ham- 
iltonian matrix  which  has  only  one  row  and  one  column.  The  matrices  representing 
the  permutations  Py  are  matrices  which  have  one  row  and  one  column.  The  matrix 
element  in  all  these  matrices  is  unity.  Therefore  the  energy  of  the  state  of  highest 
multiplicity  is  (apart  from  the  constant  in  (7.  2)  which  is  the  same  for  all  states) 
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En+1  = * IJij 

For  any  other  multiplicity  the  matrices  (P- •).  are  of  larger  size.  Let  us  consider  the 

IJ 

diagonal  energies  of  all  the  states  of  a given  multiplicity.  These  energies  are  the  di- 
agonal elements  of  (7.  2) 


Jt  - Y - J.  .(P  .) 
<x<r  Af  ij  lie 
ij 


From  the  unitary  nature  of  the  matrices  (P^)^  we  know  that 

<piiU  « 1 

Therefore  it  is  clear  from  this  fact  and  the  fact  that  J. ^ > 0 that 


"n  + 1 


We  have  therefore  proved  that  there  is  no  state  lower  in  energy  than  the  state  of  highest 
multiplicity.  One  might  ask  if  going  through  the  process  of  diagonalizing  the  Hamilton- 
ian matrix  (7.  2)  will  push  some  state  lower  than  the  state  of  highest  multiplicity.  This 
is  not  the  case  since  in  diagonalizing  the  matrix  (7.  2)  we  use  a unitary  transformation 
which  does  not  destroy  the  unitary  property  of  (P. .).  . This  unitary  property  was  all 

1 J K(T 

that  was  necessary  to  prove  the  theorem.  In  general,  in  order  to  get  some  state  other 
than  the  state  of  highest  multiplicity  lying  lowest  it  is  necessary  to  take  states  into  ac- 
count which  occupy  the  orbitals  in  a different  manner  possibly  filling  some  of  our  n or- 
bitals twice  and  leaving  others  empty. 
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3.  ON  THE  02  GROUND  ELECTRONIC  WAVE  FUNCTION 


Microwave  absorption  measurements  on  O1^  O*7  gas^  exhibit  a hyperfine 
structure.  The  observed  spectrum  agrees  with  that  predicted  on  the  basis  of  a perturb- 
ing Hamiltonian  of  the  form  bT  • ft  + cl  S * ^ is  the  electron  spin  angular  momentum 

1 7 z z 

of  the  molecule,  and  T is  the  O*  nuclear  spin  (I  = 5/2).  According  to  a ineory  by 
Frosch  and  Foley,  experimental  values  of 


<*Vr3>  ■ IC  C°;rV  ~ 1 ) P<r)  d? 


(8.  1) 


and  p(0)  can  be  deduced  from  those  of  the  parameters  b and  c.  The  angle  between  the 
internuclear  axis  (the  z-axis)  and  the  radius  vector  r from  the  O17  nucleus  to  an  elec- 
tron  is  denoted  by  x.  and  P2(x)  = (l/2)(3x  - 1).  <P2/r  > is  the  value  of  P^/r  aver- 

aged with  respect  to  the  number  density  p(r)  of  electrons  of  unpaired  spin  of  the  mole- 
cule, and  p(0)  is  the  value  of  p at  the  nucleus.  In  the  following  we  compare  the 
experimental  values  - 1.  32  and  0.  187  (in  atomic  units)  of  <P,/r3>  and  p(0),  respectively, 

obtained  by  Miller  and  Townes'1,  a'  with  values  determined  by  Meckler's  O,  normal 

(3)  L 

electronic  wave  function.  ' ' 

<P,/r3>may  be  defined  more  precisely  as 


<Pz/r3>  = J i*  (Z-z(l°l-  J-)2sjz  $ dr 


(8-2) 


r .“ 

J 


where  Sj  is  the  spin  angular  momentum  of  the  electron  measure  in  units  of  )(.  The 
summation  is  over  the  N electrons  of  the  molecule,  and  the  integration  is  ever  the  co- 
ordinates of  all  electrons  (and  is  to  include  summation  over  spin  coordinates).  J is  a 
normalized  electronic  wave  function,  the  solution  of  SchrSdinger's  equation  for  C>2  with 
a spin  independent  Hamiltonian  and  with  fixed  nuclei.  J is  an  eigenfunction  of  ifZ  and 

S with  eigenvalues  S(S  + 1)  and  Mc  = S;  S = £.s\.  Upon  integration  over  the  coordinates 
* ° J J 

of  all  electrons  but  one,  and  summation  over  all  spin  coordinates,  (8.  2)  reduces  to  (8.  1) 

in  which  p(r ) is  defined  by 


dX„ 


(8.  3) 


p(r)  = N JV(X,  X2 XN)  2sz  $(X.  X£ XN)  d«r  dX£  . . . 

and  satisfies  Jp(r)  d?  = Ji*  2SZ  J dr  = 2S.  In  case  J is  expressed  as  a linear  com- 
bination of  determinants  of  orthonormal  orbitals,  p(r ) is  a quadratic  form  in  the  orbitals, 
the  cross-terms  (overlap  terms)  arising  from  products  of  determinants  differing  by  one 
orbital. 


-39- 


(ON  THE  Oz  GKOUND  ELECTRONIC  WAVE  FUNCTION) 


Quite  generally  p(r  ) can  he  employed  to  find  the  expectation  value  of  F: 

I'i*  F if  dr  - J f(r)  p(r)  dr,  where  F is  an  operator  of  the  type  F = f(r  ) 2s  . 

(4 \ J J J z 

E.  g.  . in  the  theory  of  neutron-electron  magnetic  scattering'  ' such  an  operator  occurs 

withf(r)  = — e'*1'  ‘ r where  K is  the  neutron  scattering  vector. 

2 3 - (3 ) 

The  normal  O,  electronic  level  is  a ) Meckler'  ' approximates  $ for 

‘ g 3 

this  level  by  a linear  combination  of  nine  symmetry  states  of  T symmetry,  each  a 

‘'g 

linear  combination  of  determinants  of  orthonormal  orbitals,  these  orbitals  being  sym- 
metry orbitals  (LCAO  molecular  orbitals)  from  one  configuration.  We  write  Meckler’s 
if  as 


rM 


= Y.  c ♦ 

p = a 


(8.4) 


where,  in  terms  of  the  determinants  specified  in  Table  8-1,  the  symmetry  states  are 
given  by 


<t>  = i(A  + B - C - D) 

« im 

4>b  = i(B  + C - A - D) 

<*>  = 2'  3/2{E  + G-  F -H  + J + K-  I-L)  (8.5) 

4>h  = 2'  3/Z(F  + G-  E-  H + I + K-  J-L) 

4>i  = 2'  3/2  (E  + F - G - H + I + J - K - L) 

and  the  one-determinant  symmetry  states  4>e>  and  fy.  The  coefficients  c^  as  a 

function  of  R are  given  by  Meckler!3^  The  result  of  inserting  (8.4)  into  (8.  3)  is,  in 
terms  of  Meckler:s  symmetry  orbitals  (see  Table  8-1), 

p(r)  = Cj|*0|2  ♦ C2|x0|2  ♦ C3|xt!2  ♦ C4,*t,2  (8.6) 


where 

Cj  = a + P+  y-6 
C^  = o + p - y+  6 

c3  = lca  + cb|2  + 2lccl2  + 2|ce!2  + q - p + Y + 6 
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Table  8-1 

Determinants  of  the  symmetry  states  * 


2P*u 

Zpw 
H g 

Zp  TT 

g 

2P"u 

2P"u 

♦o 

*o 

x+ 

X. 

A 


+ - + - + 


+ 


♦ 


B 


+ - + - + 


+ 


+ 


C 


+ - + - 


+ 


+ 


+ 


D 


+ - + - + 


+ 


+ 


E 


+ 


+ 


+ - 


♦ - 


F * 


+ 


+ - + - ♦ 


G 


+ 


+ 


+ - + - + 


H 


+ 


+ 


+ - + - ♦ 


I + + 


+ 


J 


+ 


♦ 


+ - + 


+ - 


K 


+ 


+ - + 


+ 


+ - 


L 


+ 


+ - + 


+ 


+ - 


♦c 

♦d 

*e 

♦f 


+ - 


+ - 
+ - 


+ + + - 

+ - + - + 

+ + + - 

♦ - + - + 


+ - 
+ 

+ - 
+ 


*The  normalized  determinants  A,  . . . , with  Mg  = I are  defined  by  the  symmetry  or- 
bitals on  the  principal  diagonal;  these  consist  of  filled  ls<r  and  Zs<r  orbitals  accounting 
for  eight  of  the  sixteen  electrons,  the  remaining  orbitals  being  specified  by  the  above 
table,  a + or  - entry  indicating  that  the  space  orbital  is  occupied  by  an  electron  with 
+ or  - spin.  The  space  orbitals  are  labeled  both  in  standard  notation  and  in  Meckler's 


3-, 

notation.  Meckler  lists  the  nine  2. 
eigenvalue  Mg  = 0;  these  yield  the  $ 


g 


symmetry  states  as  eigenfunctions  of  Sz  with 

when  operated  on  by  (S  + iS  )/  %/T. 

x y 
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C^  = |ca  - c/  + 2|c^|^  ♦ 2|Cjl2  - q + P+  'Y  + 6 


= W ' '■>  * Ci'2 

p = ?'•  cg*ch* 

ci|2 

= K 4 ch  - s'2 

6 = i'cg  * ch  * 

M2 

No  overlap  terms  appear  in  p since  all  determinants  involved  differ  by  at  least  two  or- 
bitals. This  follows  from  the  fact  that  no  two  space  orbitals  available  for  occupation 
outside  the  lsor  and  2s<r  filled  shells  are  of  the  same  symmetry.  Values  of  the  C's  and 
their  component  parts  are  given  in  Table  8-2.  At  normal  internuclear  distance  R = 2.  28 


Table  8-2 

Values  of  the  C's  of  Eq.  (8.  6)  and  their  component  parts  for  Meckler'r 


O^  normal 


Y electronic  wave  function.  Note  that 
R 


U ck  ■ 2- °°- 


R 

2.  0 

2.  5 

4.  0 

00 

C1 

- 0. 0018 

- 0.  0036 

0.4424 

1/2 

C2 

0. 0024 

0. 0066 

0.  4425 

1/2 

C3 

1.  9632 

1.  9183 

0.  5960 

1/2 

C4 

0. 0361 

0.  0788 

0.  5181 

1/2 

a 

0. 0003 

0.  0013 

0. 2226 

1/4 

p 

0. 00002 

0.  0002 

0. 2202 

1/4 

Y 

0. 0060 

0.  01 13 

0. 00004 

0 

6 

0. 0081 

0. 0164 

0. 00008 

0 

21c  I2 
c 

1. 9400 

1. 8694 

0.  3804 

1/4 

2lcdl2 

0. 0216 

0.  0479 

0. 3208 

1/4 

2|cel2 

0. 0084 

0. 0183 

0. 2129 

1/4 

-i‘%12 

0. 0007 

0.  0041 

0. 1998 

1/4 

2 

|ca  * cb' 

0.  00000 

0. 00004 

0. 00001 

0 

|ca  + cb'2 

0.  0005 

0. 0018 

0.  0001 

0 

-42- 


(ON  THE  02  GROUND  ELECTRONIC  WAVE  FUNCTION) 

the  symmetry  state  $ , which  corresponds  to  the  molecular  orbital  solution  for  the 

c 2 

normal  O,  molecule,  predominates  (l  - |c  I « 1).  The  p for  4 , which  involves  only 

C 6 C 

2pw  orbitals  x+>  is  not  much  different  from  the  p for  jpM. 

8-3  M 

The  value  of  <P2/r  > calculated  from  p of  at  equilibrium  internuclear 

separation  is  only  about  27  percent  of  the  experimental  value  - 1.  32.  When  Meckler's 

orbitals  in  p are  replaced  by  Hartree-Fock  orbitals,  the  value  of<P,/r3>  is  about  90 

3 ^ 

percent  of  the  experimental  value.  These  values  of  <P2/r  > are  discussed  in  more 

detail  in  the  next  paragraph.  It  seems  unlikely  that  if  Hartree-Fock  or  other  reasonable 

types  of  atomic  orbitals  were  used  in  Meckler's  calculation  in  place  of  his  Gaussian- 

type  orbitals,  the  C's  of  (8.6)  would  be  significantly  different.  The  value  of  <P2/r3> 

will  then  depend  largely  on  the  types  of  orbitals.  One  expects  that  calculated  values  of 

the  molecular  energy  and  <P2/r^>  as  a function  of  the  atomic  orbitals  employed  are 

not  very  well  correlated,  since  <P,/r3>  , but  not  the  energy,  is  sensitive  to  values  of 

C 17 

the  orbitals  primarily  in  the  vicinity  of  the  O nucleus. 

Let  us  consider  first  <P,/r*>  calculated  for  $ . We  then  have,  in  terms 

17  ^ j £ 

of  the  atomic  orbitals  p+  about  the  O1 ' nucleus  and  p+*  about  the  O nucleus 

P 2 

<P2/r3>  = zj  -i  |X+ J2  dr”  * (2  • !L)  A [l  - c^R)  + t3(R)]  (8.7) 

where 

A = /<p2/r3>  tP+!2 

c2(R)  = zf  (P2/r3)  pt*  p+'dr/A  (8.7') 

«3(r)  = /(P2/r3))p+'|2  dr/A 


and  K(R)  is  a normalization  constant.  A,  which  is  - <P,/r  > for  an  isolated  O atom 

4 3 & 

in  its  ground  level  2p  P,  is  independent  of  R;  it  is  negative,  since  the  angular  part  of 
the  integral  is  negative.  For  Meckler's  Gaussian  type  p orbitals,  A is  - . 304  and,  if 
we  ignore  the  other  factors  in  (8.7),  yields  a theoretical  value  of  <P2/r3>  only  23  per- 
cent of  the  experimental  value.  The  theoretical  value  is  increased  somewhat  by  the 
remaining  factors  in  (8.  7).  The  contributions  *2  and  i3  involving  the  orbital  are 
small,  and  should  decrease  as  R increases;  in  fact,  for  Meckler's  orbitals  and  R = 2, 

«2  = - 0.  03  and  l£3/t2|  « 0.  04,  so  that  the  factor  (1  - c2  + c3)  increases*  the  theoretical 

Actually  t?  is  the  algebraic  sum  of  three  terms  whose  magnitudes  for  R = 2 are  of  the 
order  ten  times  lc2l,  and  it  is  not  obvious  that  It^l  decreases  monotonically  with  in- 
creasing R. 
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U 


value  by  about  3 percent.  K for  Meckler's  orbitals  has  the  values  1.  25,  1. 09.  and 
1.  00  at  R = 2.  0.  2.  5 and  4.  0,  respectively,  so  that  the  normalization  of  \ increases 
the  theoretical  <P,/r  > by  about  15  percent  at  normal  internuclear  separation.  The 
net  result  is  ths-.t  <P^/r  > for  the  molecular  orbital  ground  state  of  O^  computed  us- 
ing Meckler's  Gaussian  type  orbitals  is  still  or.ly  27  percent  of  the  experimental  value. 

A for  a Hartree-Fock  O orbital^a^  is  71  percent  of  the  experimental  . If, 

however,  one  includes  the  other  factors  of  (8.  7),  assuming  these  are  unchanged  in  value, 
the  discrepancy  can  be  reduced  to  about  10  - 15  percent,  perhaps  even  less  when  these 
factors  are  properly  computed.  The  value  of  <P2/r3>  for  is  in  slightly  poorer 
agreement  with  experiment  than  that  for  Now  p involves  integrals  like  (8.  7)  with 

the  2pw  , 2pa  , and  2p<r  orbitals.  Since  the  latter  integrals  contribute  to  <P?/r^> 

U g u c 

with  small  coefficients,  they  need  not  be  investigated  with  great  precision.  For  2pw 

2 2 u 
the  integral  differs  from  (8.  7)  in  that  K is  replaced  by  Meckler's  L < 1 and  -•>  - c^, 

and  so  is  smaller  in  magnitude  than  (8.  7).  For  2pir  orbitals  the  sign  of  the  integral  is 

reversed,  A being  replaced  by  - 2A.  It  is  clear  that  <P^/r3>^^  does  not  differ  appre- 

ably  from  <P,/r*>^  • 


ci 


Let  us  next  compare  the  experimental  and  calculated  values  of  p(0).  For 

p( 0)  is  ve.-y  much  less  than  the  experimental  value.  This  is  not  surprising,  since  only 

17 

symmetry  orbitals  containing  s orbitals  have  an  appreciable  density  at  the  O nucleus, 
and  involves  s orbitals  outside  the  filled  Iso-  and  2s<r  shells  only  in  the  2p«r  orbitals, 
and  then  only  in  virtue  of  the  orthogonalization  of  the  2pir  to  the  ls«r  and  2s*  orbitals. 
Table  8-3  shows  values  of  p(0)  for  some  atomic  orbitals.  Values  of  +ls(0)|^  and 

Table  8-3 

Values  of  p(0)  for  some  atomic  orbitals 


° (5)  °2  (3) 

Hydrogen-Like  Hartree-Fock'  ' Meckler'  ' 


Is 

2s 


0.  318  Z 


0.  318  (Z /Z)* 


148.  2 
7.  604 


45.  4 
0.  363 


|4i2S(0)l  for  Meckler's  orbitals  are  smaller  by  factors  0.  306  and  0.048,  respectively, 
than  the  corresponding  values  for  the  Hartree-Fock  O orbitals  of  Hartree,  Hartree 
and  Swirles.  ^ Even  when  the  latter  replace  Meckler's  atomic  orbitals  in  his  symmetry 
orbitals,  p(0)  is  less  than  0.001  for  ip^.  smaller  by  a factor  of  about  100  than  the  ex- 
perimental value  0.  187.  (We  have  neglected  the  reorthogonalization  and  renormaliza- 
tion that  should  be  performed  for  these  substituted  atomic  orbitals.  ) 

Apparently  the  normal  electronic  state  contains  small  admixtures  of  other 
symmetry  states  than  those  considered  by  Meckler,  and  these  additional  symmetry 
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states  have  large  effects  on  p(0).  The  situation  can  be  clarified  qualitatively  to  some 
degree.  A rough  criterion  for  the  magnitude  of  the  amplitude  c^  of  a symmetry  state 
is  the  sum  of  the  energies  of  the  component  orbitals,  the  higher  the  energy,  the  smaller 
I c | - Only  the  <r  oroitals  of  a symmetry  state  contribute  appreciably  to  the  density  of 

**  17 

unpaired  spin  in  the  vicinity  of  the  O nucleus,  since  only  these  contain  s orbitals. 
Thus,  only  symmetry  states  with  unpaired  a orbitals  can  contribute  effectively  to  p(0). 
We  can  get  some  idea  of  the  relative  effects  on  p(0)  of  various  nsa  orbitals  by  noting 
that  for  hydrogen-like  atomic  orbitals  4^.  I *^no  (0)/*4»  1 0 (°)I ^ = n~  ^ ' 

On  the  basis  of  these  ideas  we  may  expect  that  only  symmetry  states  with  low 
energy  containing  unpaired  nstr  orbitals  with  small  n give  important  contributions  to 
p(0).  Examples  of  such  are  symmetry  states  which  differ  from  one  of  the  in  that  a 
ls<r  or  2str  orbital  is  replaced  by  an  nl<y  of  low  n £ 2 and  If  0,  or  in  that  a 2pa  is  re- 
placed by  an  nso-  with  small  n > 3.  (Note  that  if  the  excited  symmetry  state  is  to  have 
the  same  symmetry  as  a ^ and  involve  determinants  which  differ  by  only  one  orbital 
from  thos  of  the  the  excited  orbital  must  have  the  same  symmetry  as  the  unexcited. 
Also,  a cross  term  in  p(r)  vanishes  at  r = 0,  unless  the  orbitals  which  overlap  are 
both  s a ) 


Approximate  wave  functions  in  molecular  theory  are  commonly  determined 
by  a variational  procedure  which  makes  the  energy  stationary.  The  result  is  that  al- 
though the  energy  is  insensitive  to  small  variations  in  the  wave  function,  other  quanti- 
ties, such  as  <P2/r3>  and  p(0),  may  be  relatively  quite  sensitive.  One  should  then 
not  expect  amounts  of  mixing  of  symmetry  states  which  yield  energies  in  satisfactory 
agreement  with  experiment  to  necessarily  yield  other  molecular  quantities  in  such  satis- 
factory agreement;  p(0)  is  an  example  of  such  a quantity.  On  the  same  basis,  types  of 
atomic  orbitals  which  lead  to  reasonable  energies  may  lead  to  poor  values  of  other  mo- 
lecular quantities:  both  <P2/r3>  and  p(0)  are  examples  of  such  quantities. 


References 

1.  S.  L.  Miller  and  C.  H.  Townes,  Phys.  Rev.  9£,  537  ( 1 953) . 

la.  S.  L.  Miller,  C.  H.  Townes  and  M.  Kotani,  Phys.  Rev.  90,  542  ( 1 953). 

2.  R.  A.  Frosch  and  H.  M.  Foley,  Phys.  Rev.  88,  1337  ( 1 952) . 

3.  A.  Meckler,  to  be  published  in  J.  Chem.  Phys;  see  also  A.  Meckler,  Thesis, 

Department  of  Physics,  M.I.T.  ( 1 952)  and  previous  issues  of  this  Report. 

4.  O.  Halpern  and  M.  H.  Johnson,  Phys.  Rev.  ^5,  898  (1939). 

5.  D.  R.  Hartree,  W.  Hartree  and  B.  Swirles,  Trans  Roy.  Soc.  (London)  A238,  229 

(1939).  


W.  H.  Kleiner 


-45- 


9 THE  ETHANE  MOLECULE 


As  indicated  in  the  last  Report,  the  severe  restriction  on  the  amount  of  con- 
figuration interaction  to  be  considered  necessitates  an  optimization  of  the  molecular 
orbitals,  a best  choice  of  linear  combinations  of  atomic  orbitals.  This  problem  is  faced 
also  by  Dr.  Kaplan  in  his  treatment  of  the  NK^  molecule  where  even  though  quite  a bit 
of  configuration  interaction  is  handled,  the  amount  is  still  not  enough  to  permit  complete 
equivalence  of  all  linear  combinations  of  atomic  orbitals.  Roothaan  ^ has  given  a pro- 
cedure for  determining  best  LCAO's  for  the  case  of  a closed  shell  state,  a single  de- 
terminant many-electron  wave  function.  His  equations  are  the  algebraic  equivalent  of 
the  Hartree-Fock  differential  equations  and  are  similarly  characterized  by  the  dis- 
couraging prospect  of  a self-consistent  calculation.  However,  there  is  Whirlwind  -- 
the  high  speed  computer  now  of  large  capacity  --  and  Roothaan's  procedure  can  be  mech- 
anized and  performed  by  the  machine  so  that  it  becomes  a fast,  repeatable  routine.  We 
wish  to  present  •->  flow  chart  of  the  machine  program  now  being  developed  but  first  we 
will  re-express  Roothaan's  equations  in  a slightly  modified  form: 

To  begin  with,  there  are  the  one -electron  functions  which  are  to  be  combined 
in  the  best  possible  ways.  Let  them  be  orthonormalized.  Here  we  differ  with  Roothaan 
who  did  not  require  orthogonality.  The  lack  of  orthogonality  leads  to  a more  compli- 
cated type  of  determinantal  equation  to  be  solved  and  we  see  no  point  in  adding  this  diffi- 
culty to  the  non-escapable  one  of  self-consistency.  We  begin  with  a well-considered 
first  guess  to  the  molecula.'  orbitals,  an  orthonormalized  set.  A general  linear  com- 
bination of  these  functions  is  still  a general  linear  combination  of  the  original  atomic 
functions. 

This  original  set  of  functions  is  denoted  as  v^  and  the  number  in  the  set  by  N. 

If  there  are  2n  electrons  in  the  system,  we  are  to  find  n linear  combinations  of  the  v^'s 


♦i 


y c.  v 

^ ip  p 


such  that  an  antisymmetrized  many-electron  wave  function  formed  with  each  doubly 
occupied  will  be  associated  with  a minimum  value  of  the  energy.  The  expression  for 
the  energy  is 


E ■ 2 1 vH,„  ♦ i ‘VvGi>w 

where 

* /v(,)  [•  vi2  * v(1)]  v,(1)d3xi 
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= 2(pvlX(r)  - (pcrIXv) 

(pvlXe)  = f v *(l)  v (1)  ~L-  v *(2)  v (2)  d3x,  d3x? 

J K v r j ^ K a 1 c 

p = y c*  c. 

r\LV  L.  iji  IV 

i 

or  in  matrix  notation 

p = C+C  (+  = Hermitian  adjoint) 

and  the  orthonormalization  of  the  ^'s  is  implied  by 

CC+  = 1 

C is  an  n x N rectangular  matrix,  p is  a square  N x N.  The  n rows  of  C, 
which  are  the  complex  conjugates  of  the  n columns  of  C,  are  the  n lowest  eigenvectors 
of  the  matrix 


KX.  * HXr  * 


I *■ 


pv 


G*v 

|*V  X* 


That  is 


KC+  = AC+  A diagonal 


We  are  to  pick  a p,  form  K,  diagonalize  it,  selects  its  n lowest  eigenvectors,  form 
p,  re-cycle.  The  thing  to  realize  is  that  H and  G do  not  change  from  cycle  to  cycle. 
They  are  the  input  to  the  machine,  the  data  to  be  used  over  and  over  again.  The  flow 
chart  is  illustrated  on  p.  48. 

The  machine  is  to  do  matrix-vector  multiplication  (symbolized  as  pG),  matrix 
addition,  matrix  diagonalization,  magnitude  selection  and  comparison,  and  many  cycles. 
Even  with  the  type  of  factorization  by  symmetry  discussed  by  Roothaan  there  is  a lot  to 
be  stored  in  the  machine,  but  the  addition  of  more  electrostatic  storage  tubes  and  a 
magnetic  drum  to  Whirlwind  has  made  for  the  pleasant  atmosphere  of  a buyer's  market 
in  storage  space. 

The  matrix  diagonalization  procedure,  which  is  just  a sub-routine  in  the  whole 
self-consistent  process,  has  been  programmed  and  will  exist  as  a separate  routine  in 
the  Whirlwind  library.  The  present  program  can  handle  up  to  a 31  x 31  real  symmetric 
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matrix  and  a 15  x 15  complex  Hermitian  matrix.  The  rest  of  the  L.C.A.O.  determina- 
tion is  now  being  programmed  and  will  be  first  applied  to  the  NH^  molecule. 

Alvin  Meckler 
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10.  LIMITED  CONFIGURATION  INTERACTION  TREATMENT  OF  THE  NH3  MOLECULE 

The  calculation  of  matrix  elements  indicated  in  the  last  Progress  Report  is 
still  underway.  The  matrix  elements  have  been  expressed  in  terms  of  "directed"  mo- 
lecular orbitals.  Due  to  the  symmetry  properties  of  these  one-electron  functions, 
there  are  only  thirty-five  independent  two -electron  integrals  and  five  independent  one- 
electron  integrals.  Each  of  these  integrals,  however,  is  a combination  of  many  similar 
type  integrals  involving  atomic  wave  functions.  At  present  the  component  atomic  in- 
tegrals are  being  expressed  in  a form  necessary  for  the  evaluation  of  their  radial  sub- 
integrals by  I.  B.  M.  machine  methods. 


H.  Kaplan 
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11.  THE  WATER  MOLECULE 


The  tabulation  of  the  an's  mentioned  in  the  last  Quarterly  Progress  Report 
is  nearly  complete.  We  now  have  the  following  ran's; 

r = 0(.  1)  1.  2(.  2)  4.  0(.  5)  7{l)  12 

a * 1.4,  1.5,  1.6,  1.8,  2.0,  2.2,  2.6,  3.0,  4.  0 and  5. 0 

k = 1,  2 

n = 0,  1,  2,  3,  4,  5 (for  k = 1) 

n = 0,  1,2,  3 (for  k = 2). 


The  interelectronic  interaction  integrals  involves  terms  of  the  form 


(11.  1) 


These  can  be  written  as 


(11.2) 


We  shall  call  the  quantity  in  square  brackets  in  (11.2)  the  inner  potential. 

We  have  written  and  tested  a program  for  computing  the  inner  potential, 
where  is  the  product  of  two  radial  probability  amplitudes,  for  the  I.  B.  M.  Card 
Programmed  Calculation.  Our  mesh  is 


r « 0(.  0125)  . 05(.  025)  . 2(.  05)  . 6(.  1)  1.  2(.  2)  4(.  4)  10  . 


The  program  takes  three  minutes  of  machine  time  per  integration.  The  results  of  this 
indefinite  integral  are  printed  and  punched  on  cards  ready  for  use  in  definite  integrals. 
This  program  is  now  in  use. 


G.  F.  Koster 
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12.  CONFIGURATION  INTERACTION  APPLIED  TO  THE  HYDROGEN  MOLECULE 


Calculation  of  the  several  matrix  elements  discussed  in  the  previous  Progress 
Report  is  continuing.  The  odd  one-electron  function  obtained  from  the  Hartree  equa- 
tion containing  an  "average  field"  due  to  an  electron  in  the  even  orbital  has  been  deter- 
mined. The  kinetic  energy  and  nuclear  framework  integrals  appearing  in  the  main  di- 
agonal element  in  which  the  electrdns'  arein  the- odd  orbital  -have -bees. computed  aAd  . 
there  remains  only  the  electron -electron  interaction  integrals  to  be  found. 

E.  Cailen 
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13.  CONFIGURATION  INTERACTION  FOR  THE  FLUORINE  MOLECULE 


(1) 


In  order  to  shorten  the  computation  in  the  problem  mentioned  in  the  last 
Progress  Report,  it  has  been  decided  to  use  Slater-type  atomic  orbitals^  instead  of 
Hartree-Fock  atomic  orbitals.  Work  on  the  problem  is  proceeding  along  these  lines. 
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14.  NUCLEAR  ELECTRIC  QUADRUPOLE  INTERACTION  IN  THE  KC1  MOLECULE 


The  introduction  of  a polarization  effect  by  appropriate  distortion  of  atomic 
orbitals  has  been  studied  along  lines  indicated  in  the  previous  Progress  Report.  ^ A 
special  case  of  the  proposed  form: 


4- 


= Y? 


) 


(14.  1) 


has  been  investigated  quantitatively.  In  particular,  we  have  chosen  n = 1,  i = 0,  m = 
0.  /'  = 1.  m'  = 0.  and  in  addition  we  have  provided  for  normalization  and  a "mixing 
coefficient",  C.  Thus  our  wave  funciion  is: 


4*  * Ne  r(l  + Crv  cos  0) 


(14.2) 


where  N = normalization  constant  and  C,  v * variation  parameters.  (Unless  otherwise 

-8 

stated,  atomic  units,  a.  u.  , are  used  throughout;  a = unit  of  length  = .529  x 10  cm, 
2 o 

e /2aQ  * unit  of  energy  = 13.  53  e.  v. ) 

Using  (14.  2)  we  have  considered  two  types  of  perturbing  fields:  (l)  a weak 

uniform  field  and  (2)  the  field  of  a point  charge  a distance  R = 1.  4 a.  u.  away  from  the 

center  of  the  unperturbed  function. 

Case  (l) 


z 


The  Hamiltonian  is: 


H « HQ(r ) ♦ H'(r , <f  ) = - v2  - £ + / r cos  8 


and  H'(r,^  ) is  considered  to  be  small  compared  to  Hq.  We  form: 


E 


E 

o 


+ . . . 


(14.  3) 


Our  problem  is  really  that  of  the  second  order  Stark  effect  for  the  normal  hydrogen 
atom  --  which  has  been  solved  exactly^  yielding  a polarizability:  a * 2.  25  a.  u.  (4.  50  aQ^ 
in  any  other  system  of  units).  Using  (14.  2)  and  varying  C and  v we  obtain  a minimum 
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of  (14.  3)  for  C = - . 724,  v = 1.  50,  giving  a = 2.  25  a.u.  --  the  correct  value.  For 
comparison  we  have  used: 


*1  ■ »l  k * Cl*2pr] 


r cos  -r/2 
e ' 

4 ^/Zm 


(14.4) 


Minimizing  (14.  3)  with  respect  to  Cj  we  obtain  a c 1.  48  a.  u.  (2.  96  aQ^)  — a value  in 
error  by  35  percent.  The  resulting  expressions  for  and  *i*  are  given  below  for  Case 
(l).  (We  have  separately  normalized  the  angular  dependent  term  of  (14. 2)  so  that  the 
coefficient  of  the  distortion  terms  are  directly  comparable): 


+1=N1 

[+.»  - -”5'+2p,] 

ni  1 

1 - .495*2 

r 

, 1.  50  Q -r  i 

+ * N I 

,-r  - 1.4 •>*{*  221?' 

N = . 

564  - . 369  & 

L 

It  should  be  noted  that  for  Case  (l)  (14.  2)  is  just  of  the  form: 


= +°(l  ♦ AH'  ♦ . . . ) 

which  was  suggested  by  Lennard-Jones^  on  the  basis  of  an  approximate  treatment  of 
second  order  perturbation  theory.  Also  Slater  and  Kirkwood'  ' have  treated  the  normal 
hydrogen  atom  in  a weak  uniform  field  using  a function  of  the  form: 

* « +°(l  + H'R(r)) 

They  obtained  a differential  equation  for  R(r)  ard  their  resulting  wave  function  was: 

4»  B N [e'r  - \i  |(r  ♦ ^)  cos  9 e'rj]  N - .564  - . 282*  2 

and  a * 2.  25  (4.  50  a^^)  - the  exact  value- 

o 
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Case  (2) 


The  Hamiltonian  is: 

H = HQ(r)  + H'(r\  R)  « 


1 


r' 


= ^ Z pk(cos  0)  r < R 
k 

* p Y.  (p)k  *Vcos  0) r > R 

r k 


Here  we  have  a strong  field  case  with  H'(r,  R)  the  same  order  of  magnitude  as  HQ(r). 
Again  using  ( 14.  2)  and  minimizing  the  energy  with  respect  to  C and  v for  R = 1.  40  we 
find  C = . 385,  v = .65.  For  ( 14.  4)  we  find  Cj  = . 109  and  the  respective  functions  are: 


* = N je~r  + . 298  (t-15-?,?^.0-^)]  N = . 508 

L v y^oo  yj 

+ 1 - Nl[+ls  + 109  +2pz]  N1  = 989 

The  exact  solution  of  Schrodinger's  equation  for  this  potential  is  the  hydrogen  molecule 

ion  with  an  energy  at  R = 1. 40  of  - 1.  15  a.  u.  + gives  an  energy  - . 916  while  gives 

* + 

- . 809,  the  first  representing  a 20  percent  greater  value  than  the  exact  value  for  H^ 
and  the  second  a 30  percent  greater  value.  However  we  are  trying  to  solve  a polariza- 
tion problem  and  not  the  H£*  problem  and  thus  we  have  not  brought  in  the  symmetry  of 
the  positive  charges.  Our  calculation  only  shows  that  the  wave  function  of  ( 1 4.  2)  is 
more  satisfactory  than  the  case  of  (14.4)  where  we  admix  some  2pz  function  with  the  Is. 
The  method  of  Slater  and  Kirdwood  depends  for  the  separability  of  their  differential 
equation  on  linearity  of  H'  in  rectangular  coordinates  and  spherical  symmetry  of  i|*° 
and  appears  not  to  be  readily  manageable  for  complicated  cases. 

As  a further  test  of  the  function  (14.  2)  the  determinantal  wave  function  for  F” 
is  being  set  up  using  one -electron  functions  of  the  form  (14.  2)  and  determining  C and 
v by  minimizing  the  total  energy. 

References 

1.  L.  C.  Allen,  Quarterly  Progress  Report,  Solid-State  and  Molecular  Theory  Group, 
M.I.T.,  April  15,  1953,  p.  51. 


-56- 


(NUCLEAR  ELECTRIC  QUADRUPOLE  INTERACTION  IN  THE  KC1  MOLECULE) 


References  (con'd) 

2.  See  e.  g. , L.  Pauling  and  E.  B.  Wilson,  Introduction  to  Quantum  Mechanics  (McGraw- 

Hill,  1935),  p.  195. 

3.  J.  E.  Lennard-Jones,  Proc.  Roy.  Soc.  A129,  598  (1930). 

4.  J.  C.  Slater  and  J.  G.  Kirkwood,  Phys.  Rev.  .37»  682  ( 1 9 3 1 ). 

L.  C.  Allen 


15.  INTERSTITIAL  DIFFUSION  IN  THE  DIAMOND  STRUCTURE 


We  wish  to  derive  the  diffusion  coefficient  for  interstitial  diffusion  ir  the  dia- 
mond structure  by  considering  the  elementary  process  of  a transition  of  !1t  interstitial 
ion  from  its  original  site  to  a neighboring  one.  We  here  consider  only  inters. :tial  diffu- 
sion, for  the  case  of  infinite  dilution  of  the  diffusing  substance,  constant  temperature, 
and  absence  of  electric  fields. 

The  diamond  structure  can  be  viewed  as  two  interpenetrating  face-centered 
cubic  lattices,  one  having  a parallel  displacement  of  (a/4,  a/4,  a/4)  with  respect  to  the 
other,  a being  the  length  of  the  side  of  the  unit  cube  The  interstitial  positions  are 
found  to  have  the  diamond  structure  also,  displaced  by  a body -centering  translation 
from  the  original  lattice.  That  is,  if  the  carbon  atoms  are  found  at  (000),  F.  C. , and 
(a/4,  a/4,  a/4),  F.  C. , then  the  interstitial  sites  are  located  at  (a/2,  a/2,  a/2),  F.  C. , 
and  (3a/4,  3a/4,  3a/ 4),  F.  C.  The  interstitial  positions  located  with  respect  to  it  as 
as  origin  at  rTj  = (a/4,  a/4,  a/4).  n2  = (a/4,  - a/4,  - a/4),  n3  = (-  a/4,  a/4,  - a/4), 
n^  c (-  a/4,  - a/4,  a/4),  and  ihe  interstitial  position  (3a/4,  3a/4,  3a/4)  has  nearest - 
neighbor  sites  located  with  respect  to  it  by  the  negatives  of  these  displacements. 

We  assume  that  the  transition  probability  is  X.  for  a jump  from  one  site  to  a 
nearest -neighbor  site  and  O for  any  other  jump.  Then,  if  P(rj,  t)  be  the  probability 
that  site  r.  is  occupied  by  an  interstitial  ion  at  time  t,  we  have 


P(r  t)  = - 4XP(?  t)  + £ X P(?  t nr  t) 

J J i = l J 


where  the  upper  sign  holds  if  r^  is  on  the  lattice  containing  (a/ 2,  a/2,  a/2)  and  the  lower 
sign  holds  for  the  other  case.  A Taylor's  series  expansion  provides  immediately  the 
diffusion  equation 

— (?.  t)  = ^ V2  P(?,  t)  . 

at  8 

The  diffusion  constant  is  thus  D = Xa2/8. 

The  neglected  terms  in  the  equation  for  either  lattice  of  ions  are  the  third- 
derivative  terms,  but  after  averaging  over  both  lattices  the  only  terms  which  survive 
are  the  fourth  and  higher  derivatives. 

The  assumption  of  infinite  dilution  of  the  diffusing  substance  can  be  relaxed 
somewhat,  for  it  is  found  that  if  the  probability  of  occupancy  of  target  sites  is  considered 
the  same  diffusion  equation  results.  The  jumping  probability  X will  be  a function  of 
local  concentration  of  diffusing  substance,  however,  so  the  requirement  of  fairly  high 
dilution  would  still  hold. 


H.  C.  Schweinler 
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16.  A STLDY  OF  2Zp  IN  ATOMS 

The  results  of  calculations  of  the  effective  nuclear  charge  for  potential  in 
an  atom  have  been  submitted  to  the  Physical  Review  for  publication. 


A.  J.  Freeman 


17  ATOMIC  MULTIPLET  SEPARATIONS 


Atomic  multiplet  separations  have  been  computed  for  atoms  with  uncompleted 
3p  and  3d  shells.  Simple  analytic  radial  parts  of  the  one -electron  wave  functions  were 
assumed: 


- Z^r  2 *^3r 

ft3p  = C8re  + CV  e 


2 " ^4r 
R3d  = C10r  e 


R,  was  obtained  by  taking  a linear  combination  of  single  term  expressions  for  the  2p 

( i ) 

and  3p  radial  functions.  The  single  term  expressions  are  given  by  Professor  Slater  ' 
as: 


n*  - 
r 


\z  - S)/n*] 


r 


where  Z is  the  atomic  number  and  n*  and  S are  given  by  simple  rules.  Cg  and  are 
determined  from  orthogonality  conditions  with  the  2p  radial  function  and  by  normaliza- 
tion. If  one  assumes  the  radial  parts  of  the  one -electron  functions  do  not  vary  from 
multiplet  to  multiplet  for  a given  configuration,  a multiplet  separation  is  given  by  a 
sum  of  Coulomb  and  exchange  integrals.  This  sum  may  be  expressed  in  terms  of  the 
Fk(n,  I;  n\  V)  and  Gk(n',  1;  n\  V)  integrals.  ^ 

If  orthonormal  functions  of  the  form. 


R 


n I 


» 1 

k 


C.  r e 
k 


-v 


are  used,  the  integrals  may  be  written: 


W n'l')  = 2 l CsC,CuCv  kL 
stuv 


ns  + rV  nu  + nv 


(Z  + Z.,  Z + Z ) 

' S t U V7 


where: 


00  00 


k + 1 


nn.  . . n. 

L'Ua.b)=  //  r'.-'  r(b)l 

o o 

, ,,ninj9ni  + nj  ,00, 

= (-1)  — rn — rr  (a- b) 

9a  ' 9b  2 * 


,k  n 


iML  r'  j e‘br’  r*r'2 


dr  dr' 
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Multiplet 

Table  17-1 

Calculated 

Observed 

Separation 

Separation 

Separation 

Silicon  I 

3P  - *d 

10,  285.  5 cm  * 

6,  198.  66  cm 

(3p)2 

*d  - *S 

15,  428.  2 

9,  095.  43 

Phosphorus  I 

4s  - 2d 

17,  179.  0 

11,  369.  1 

(3p)3 

2d-2p 

11,  452.  6 

7,  366.  1 

Potassium  V 

4s  - 2d 

20.  914. 6 

24,  118. 5 

Op)3 

4s  - 2d 

34,  857.  8 

39,  904.  5 

3f  - *d 

6,  131.  02 

7,  069.  62 

Titanium  I , ■> 

(3d)2  (4s)2 

F - JP 

7,  099.  02 

8,  324. 8 

3F  - *G 

9,443.88 

11,  932.  79 

3F  - *D 

7,  810.  79 

10,  611.  0 

Vanadium  IV  , , 

(3d)2 

F - JP 

9,  044.  00 

12.  930.4 

3f  - !g 

12, 031. 3 

18.  040.  0 

3F  - *D 

9.  490.  51 

12.  647.  0 

Chromium  V , , 

(3d)2 

F - ■’P 

10,  988. 9 

15.  190.  0 

3F  - lG 

14,  618.  7 

21.  510.  0 

3F  - 

1 1,  170.  3 

14,  531. 

Manganese  VI 
(3d)2 

3F  - 3P 

12,  933.  9 

17,  350.  7 

3F  - !g 

17,  206.  1 

24, 706. 

3f  - 

12, 850. 0 

16,  350. 

Iron  VII  \ \ _ 

F - P 

14.  878.  8 

19,  455. 

(3d r , . 

F - ‘G 

19,  793.  5 

27,  790. 

Vanadium  I 

4 1 

4F  - 4P 

8.  363.  25 

9.  415.  24 

(3d)J  (4s) 

-61- 


(ATOMIC  MULTIPLET  SEPARATIONS) 


Table  17-1  (con'd) 


Multiplet 

Separation 

Calculated 

Separation 

Observed 

Separation 

Chromium  IV 
(3d)3 

4F 

-4P 

10,  308.  2 cm'  1 

1 3,  8 16.  cm  * 

Manganese  V 
(3d)3 

4f 

-4P 

12, 253. 2 

16, 004. 

Iron  VI 
(3d)3 

4f 

-4P 

14.  198.  1 

18, 169. 

Cobalt  VII 
(3d)3 

4f 

-4P 

16,  143.  1 

20.  271. 

Cobalt  I 
(3d)7(4s)2 

4f 

-4P 

13,420.  1 

13.  068. 98 

Nickel  I 
(3d)®(4s)2 

3f 

- 3P 

14, 684. 4 

14,  604.  17 

The  "observed  separations"  were  obtained  in  all  cases  by  averaging  the  fine  structure 
splittings  in  the  levels  concerned.  The  experimental  data  were  taken  from  Atomic 
Energy  Levels,  Charlotte  E.  Moore,  Circular  of  the  National  Bureau  of  Standards  467, 
1949. 

k 

By  definition,  r(a)  is  the  smaller,  r(b)  the  larger  of  r and  r\  The  G integral  is  ex- 
pressed as  a similar  sum.  A table  of  L functions  was  made  for  n.,  n and  k values 
from  zero  to  four,  which  is  adequate  to  evaluate  integrals  involving  electrons  up  to  the 
3d  shell. 

The  ground  state  multiplet  separations  have  been  worked  out  for  a number  of 
elements.  These  results  are  given  in  Table  17-1.  The  multiplet  separations  for  the 
(3d)2  configurations  reduce  to  linear  functions  of  atomic  number  and  are  given  by: 

3F  - 3P  * 5,  834.  8 Z4 

3F  - lG  = 7,  762.  1 Z. 

4 Z,H^ 

3 1 ^ 3 

JF  - lD  = 5,  039.  2 Z4 

3F  - lS  = 19,387.  Z4 
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There  is  good  agreement  between  the  slopes  of  the  observed  and  theoretical 

separations  as  a function  of  Z.  The  (3d)^,  *F  - *P  separation  exhibits  the  same  linear- 

7 8 

ity  in  atomic  number.  The  apparent  accuracy  in  calculating  the  (3d)  and  (3d)  multiplet 
separations  is  somewhat  accidental. 
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18.  AN  ANALYTIC  DETERMINATION  OF  ATOMIC  ENERGY  LEVELS  IN  MAGNESIUM 


Atomic  energy  levels  for  magnesium  in  various  states  of  ionization  and  ex- 
citation were  calculated  using  the  determinantal  method.  ^ The  experimental  and  cal- 
culated values  for  the  configurations;  Is2,  ls22s22pk,  ls^2s^2p^3s^,  ls22s22p^3s3p 
(*P  and  *P)  are  given  in  Table  18-1.  For  the  closed  shell  configurations  a single  de- 


terminant  was  used. 

For  the  3s3p  energy  level  a 

secular  equation  c 

orresponding  to  the 

spin-orbital  degeneracy  was  ."olved. 

Table  18-1 
Experimental 

Calculated 

Error 

Configuration 

value  (cm"  *) 

value  (cm-  *) 

percentage 

is2 

- 3.0011  x 107 

- 2.  9963  x 107 

. 15 

, 2,  2,  6 
Is  2s  2p 

- 4. 3759  x 107 

- 4.  3437  x 107 

. 74 

, 2,  2,  6,  , 
Is  2s  2p  3s3p 

3P 

- 4. 3931  x 107 

- 4. 3589  x 107 

. 78 

lP 

- 4. 3944  x 107 

- 4.  3614  x 107 

. 78 

Tne  one -electron  wave  funciions  used  in  the  determinants  were  an  orthonor- 
mal set  formed  from  linear  combinations  of  simple  analytic  shielding  constant  wave 
functions.  ^ The  normalized  radial  wave  functions  used  are  given  in  Table  18-2. 

Table  18-2 


80.  0403  e"  1 1‘  7 r 

23.  7035  e” 1 7 r - 36.  7556  e' 3>  925  r 
35. 2426  r e' 3>  925  r 

2.  57834  e”  1 *'  7 r - 4.  25902  r e'  3'  925  r + . 354533  r2  e*  * 95  r 

, nun  - 3.  925  r 2 - . 95  r 

3.  93672  re  - . 354539  r e 
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19.  THERMAL  VIBRATIONS  OF  CU-ZN  SYSTEM  CRYSTALS 


The  determination  of  the  various  branches  of  the  normal  vibration  frequency 
versus  reciprocal  lattice  vector  contour  map  and  the  associated  frequency  spectra  (num- 
ber of  normal  frequency  modes  per  unit  frequency  range)  of  the  nuclear  framework  of 
crystals  Is  of  direct  interest  in  the  thermal  diffuse  scattering  of  x-rays,  in  restrahlen 
and  Raman  spectra  interpretation,  and  in  the  calculation  of  such  thermodynamic  quan- 
tities as  specific  heats  and  coefficients  of  thermal  expansion.  There  are  three  main 
aspects  of  the  problem  as  treated  on  the  Born-von  Karman  model:  (1)  the  determina- 
tion of  the  generalized  force  constants  between  nuclei  ("dynamical  matrix"  elements 
of  Born^),  which  appear  in  the  Taylor  expansion  of  the  nuclear  potential  energy  (about 
the  equilibrium  cohesive  energy  value  at  which  all  first  derivative  or  force  terms  are 
zero)  as  coefficients  of  the  terms  quadratic  in  the  displacements  of  the  nuclei  from 
equilibrium:  (2)  the  transformation  of  the  coupled  equations  of  motion  of  the  nuclei  into 
the  normad  coordinate  form  of  linear,  homogeneous,  algebraic  equations  in  the  ampli- 
tudes of  the  normal  coordinates;  (3)  the  solution  for  the  normal  frequencies  of  the  re- 
sulting consistency  determinants  of  order  three  times  the  number  of  atoms  per  unit 
cell,  and  the  calculation  of  the  frequency  spectrum  either  from  this  frequency  plot  or 
directly  from  these  secular  determinants. 

Inherent  ir  the  method  in  the  standard  form,  of  course,  is  a neglect  of  quar- 
tic  and  higher  order  displacement  terms  in  the  potential  energy,  which  would  lead  to 
anharmonic  effects.  ^ Step  (2)  above  is  equivalent  to  the  well-known  problem  of  diag- 
onalization  of  quadratic  forms^  (here  the  simultaneous  diagonalization  of  the  kinetic 
energy  quadratic  terms  in  the  velocities  and  the  potential  energy  quadratic  terms  in  the 
coordinates),  and  if  the  unit  cell  is  simple  this  step  presents  little  formal  difficulty. 
Much  attention  has  been  focused  on  step  (3),  in  which  the  sampling  procedure  of  Black- 
man, ^ the  moment  method  of  Montroll,  ^ and  various  graphical  and  mechanical 
methods  for  obtaining  the  frequency  spectrum  have  all  been  used.  The  recent  work  of 
von  Hove^  on  the  topology  of  the  frequency  contour  map  has  provided  a basic  under- 
standing of  the  general  features  of  the  frequency  spectrum.  Quite  often  step  (1)  is  not 
handled  with  as  much  care  is  the  other  two.  Where  real  crystals  are  treated  the  force 
constants  are  usually  determined  from  experimental  values  of  the  elastic  constants. 

In  order  to  reduce  the  number  of  independent  parameters  to  the  number  of  independent 
elastic  constants  and  for  general  simplification,  a number  of  assumptions  are  often 
made:  namely,  those  of  crystal  isotropy;  centred  forces;  the  potential  being  a sum  of 
interaction  energy  terms  between  atoms  taken  only  two  at  a time;  and  the  vanishing  of 
force  constants  for  neighbors  further  away  than  second  or  third  nearest  (save  for  ionic 
cry«al.<7>).<8> 

In  the  present  treatment  of  Cu-Zn  system  crystals  particular  attention  is 
being  paid  to  the  problem  of  determining  the  force  constants  in  as  general  and  rigorous 
a way  as  possible.  The  symmetry  properties  of  the  crystals,  as  formulated  in  terms 
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of  the  well-known  theory  of  the  representations  of  space  groups,  v ' are  of  great  help  in 
simplifying  this  problem.  In  obtaining  the  remaining  independent  force  constants 
considerable  use  is  to  be  made  of  the  experimental  results  being  obtained  on  the  fre- 
quency contour  maps  of  these  crystals  in  the  work  of  the  x-ray  group  of  Professor  B. 

E.  Warren  at  M.  I.T.  cn  the  temperature  diffuse  scattering  of  x-rays/*’ 
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